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Abstract 

We consider two lattice models for strongly correlated electrons which are exactly- 
solvable in one dimension. Along with the Hubbard model and the su(2|2) spin chain, 
these are the only parity-invariant models that can be obtained from Shastry's R-matrix. 
One exhibits itinerant ferromagnetic behaviour, while for the other the electrons form 
bound pairs and at half-filling the model becomes insulating. We derive the TEA equa- 
tions for the models, analyze them at various limits, and in particular obtain zero tem- 
perature phase diagrams. Furthermore we consider extensions of the models, which 
reduce to the Essler-Korepin-Schoutens model in certain limits. 
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1 Introduction and Summary 

There exist a range of materials whose behaviour is strongly correlated and for which a 
proper understanding has yet to be established. Prominent in this list are the high-Tc 
superconductors, itinerant ferromagnets and heavy fermion systems. In this paper we 
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present models which we beheve capture some of the features of the first two, and will 
also briefly comment on a possible connection to the third. 

The complete description of a solid is a complex many body problem. Simplifications 
and exact methods are crucial for progress to be made. In particular, exactly solvable 
models provide a bedrock on which more realistic models can be constructed and ex- 
amined. The most important such model is that of free particles. It can be solved in 
any dimension for both bosons and fermions. The physics of the majority of metals 
is understood by considering the effects of interaction perturbatively on an exact free 
solution. 

Strongly correlated materials however are those which cannot be understood by re- 
ducing the complexity to a non-interacting picture. The collective behaviour of the 
constituents can drive the material into unexpected phases and interacting exactly solv- 
able models are invaluable for getting insight into such effects. The Hubbard model 
[1] is a prominent example. Indeed it has become a cornerstone of investigations into 
strongly correlated electronic behaviour. The model is integrable in one dimension, and 
describes a Mott insulator [2], which can be seen through its exact solution [5]. The 
book [1] provides an extensive review of the one- dimensional Hubbard model. 

The primary focus of this paper is to investigate two models closely related to the 
Hubbard model, we call them the A-model and the B-model. Like the Hubbard model 
they describe spin-1/2 electrons interacting on a lattice and in one dimension can be ob- 
tained from Shastry's R-matrix [5] . For this reason we refer to them as Hubbard-Shastry 
models, and by construction they are exactly solvable in one dimension. Normally an 
R-matrix depends on two spectral parameters through their difference, and is thus in- 
variant under shifts of the spectral parameters. Shastry's R-matrix is special however 
as it is of a non- difference form and so the actual value of the spectral parameters are 
important. This was noted in particular by [6l [7], who used this extra variable to con- 
struct a one-parameter extension of the Hubbard model. Their extension is hermitian 
for purely imaginary values of this parameter and is in fact also hermitian along another 
line for complex values of the parameter. The model is parity invariant however only 
for special distinct values of the parameter and these correspond, in addition to the 
Hubbard model, to the A- and B-models, and the su(2|2) spin chain. 

The A- and B-models have symmetries that are similar to those of the Hubbard 
model but the phases they exhibit are different. The A-model ground state contains spin- 
polarised electrons and thus describes itinerant ferromagnetism. For the B-model the 
electrons form bound pairs and at half-filling it becomes insulating. This is reminiscent 
of the behaviour of the supersymmetric t-J model [8], and has relevance for the modelling 
of high temperature superconductivity. 

The A- and B-models are more involved than the Hubbard model. Both have cor- 
related hopping terms, where the magnitude of the hopping parameters depend on the 
occupation of the sites by other electrons, and take into account processes such as nearest 
neighbour Coulomb interaction, spin exchange of two neighbouring electrons of opposite 
spin, and the pair hopping of two electrons from one site to a neighbouring site. The 
relative contributions of these terms are strongly constrained by the integrability of the 
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models. Like the Hubbard model they both have one free coupling constant u after one 
normalises the Hamiltonian so that hopping amplitudes are of order 1. 

The A-model does not have a free fermion limit, rather the weak coupling u — ?• limit 
is a model of graded permutations, the 5u(2|2)-spin chain [S]. In the strong coupling limit 
the Hubbard on-site interaction is the dominant term in the Hamiltonian. The model 
behaves as a ferromagnetic t-J model however, in contrast with the antiferromagnetic 
behaviour of the Hubbard model in this limit. 

The B-model does describe free spin-1/2 fermions in the weak coupling limit. As the 
coupling is increased the hopping becomes correlated and all of the interactions men- 
tioned above begin to contribute. The Hubbard interaction appears in the Hamiltonian 
but there is no range of u for which it is dominant. In the strong coupling limit all 
terms are of the same order of magnitude, and the model appears to be equivalent to 
the Essler-Korepin-Schoutens (EKS) model with coupling U = A. Furthermore, the 
model possesses a hidden su(2|2) fermionic symmetry for arbitrary coupling. 

The integrability of the models is used to study their exact solution in one dimension. 
As for the Hubbard model, this can be done for arbitrary chemical potential ^ and 
applied magnetic field B. Our formalism allows us to study the Hubbard, A- and 
B-models in parallel. We restrict our attention to the parameter region < and 
i? > (results for the other regions can be obtained using symmetries), here the electron 
number density ranges from (empty lattice) to 1 (half-filling) and the magnetisation 
is non-negative. The TEA equations of the models are derived and we analyze them for 
various limits. Particular focus is paid to obtaining the zero temperature phase diagram 
in the thermodynamic limit, as this provides the most insight into their properties. 
Diagonalising the models using the Bethe ansatz technique, taking the thermodynamic 
limit, and going to zero temperature, the ground states are found to be composed 

• y-particles: spin-up electrons, 

• w-particles: each change the spin of the state by -1, 

• l|fw-strings: each is a bound state of one spin- up and one spin-down electron. 

It is well established |1] that for the Hubbard model there are no l|ww-strings in the 
ground state. In contrast for the B-model there are no w-particles, and with no applied 
magnetic field there are only l|f w-strings. The ground state of the A-model is composed 
only of y-particles. The phase diagrams for the A- and B-models are presented in Figure 
[ij The phases identified are: I. empty band, II. partially filled and spin polarised band, 
III. half-filled and spin-polarised band, IV. partially filled and partially spin polarised 
band, V. half-filled and partially spin polarised band. 

The ground state of the A-model is spin polarised as it contains only spin-up electrons 
for B > Q. The electron number density increases for increasing chemical potential, the 
lines separating phases I and II and phases II and HI depend on the coupling constant 

^ In the conventions of [3] a w-particle is a A-string of length 1, a l|tiw-string is a fc-A string of 
length 1, and y-particles could have been called fc-particles. 
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Figure 1: Zero temperature phase diagrams in the ^S-pIane for the A- and B- and extended 
B-models. The phases identified are: I. empty band, II. partially filled and spin polarised band, 
III. half-filled and spin-polarised band, IV. partially filled and partially spin polarised band, 
V. half-filled and partially spin polarised band. The dependence of /xo,//i,//2 on u is given in 
the text, /Lie = —2 -|- 2 log 2, and g is the coupling constant of the extended B-model. 



respectively as yUi = -2 - 2^1 + - B and yU2 = 2 - 2^1 + - B. For i? < 
the phase diagram is mirrored except that here the ground state is polarised with spin- 
down electrons. Thus we see that the A-model, which describes itinerant electrons on a 
one-dimensional lattice, is ferromagnetic. 

The phase diagram of the B-model is richer. In the weak coupling limit the point 



/io, whose dependence on u is given by eq.(6.6), approaches 0, phase V shrinks to a 
portion of the 5-axis, and the phase diagram reduces to that of free electrons just as 
the Hamiltonian does. As the coupling constant is increased to its strong coupling limit 
the point /zq decreases monotonically to the value /Zc = — 2 -|- 2 log 2 ^ —0.6137. Let us 
first examine the phase diagram as increases along the line 5 = 0. Here the ground 
state has zero magnetisation and contains only bound pairs of spin-up and spin-down 
electrons (the l|f ly-strings). For yU < —2 the cost of having particles in the ground state 
is too high and the band is empty. The ground state then fills up with the l|f ti'-strings 
as 11 is increased from —2 to /^o, and so is a liquid of singlet pairs. When jj, reaches 
/io the ground state is half-filled and remains half-filled for a range of /x. Thus a gap 
opens for charge excitations of the ground state and phase V is insulating. Consider 
now i? > 0. As i? is increased spin-up electrons enter the ground state and the model 
becomes magnetised. This increases the value of /i that separates phases IV and V and 
eventually spin-polarises the model, breaking up all of the bound pairs. 

It is possible to extend both the A- and B-models by adding the Hubbard interaction 
term to their Hamiltonians. This decouples the Hubbard interaction from the other 
terms and in particular for the B-model it allows strong onsite Coulomb repulsion of 
electrons to be introduced. When the extension dominates the B-model reduces to the 
t-J model with coupling J interpolating between and 2 as u changes from to oo. 

In general the extension destroys the integrability of the models but it happens that 
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there are special values of the coupling constant u for which it does not. If u = the 
extended B-model is just the Hubbard model. In addition the EKS model for arbitrary U 
can be obtained from the B-model when the coupling constant goes to infinity. Similarly 
the extended A-model reduces to the EKS model up to the sign of the Hamiltonian after 
taking u to 0. Our formalism can be naturally adapted to the EKS model, and we derive 
the simplified TEA equations and analyze the zero temperature phase diagram including 
the dependence on applied magnetic field, see Figure [l] 

The plan of the paper is as follows. In section 2 the models are introduced on the 
level of Hamiltonians and their symmetries are presented. The Bethe equations for the 
Hubbard-Shastry models are discussed in section 3, and the parity invariant models are 
identified. Section 4 outlines our formalism for the TEA analysis of the thermodynamic 
properties of the models, in particular the string hypothesis is formulated, and the TEA 
equations and free energy are presented. In section 5 we investigate the TEA equations 
in various limits, namely the weak coupling, strong coupling, strong magnetic field, 
large chemical potential and infinite temperature limits. The zero temperature limit is 
discussed in detail in section 6 where the phase diagrams for the A- and E-models are 
obtained. In section 7 the integrable extensions of the A- and E-models, equivalent to 
the EKS model, are studied. In Conclusion we discuss future directions. In the appendix 
A.l our conventions, definitions and notations are summarised. The derivation of the 
free energy and the canonical and simplified TEA equations is outlined in appendix A. 2. 
The relation between repulsive and attractive models is discussed in appendix A. 3 on 
the level of the TEA equations. The algebraic limit of the TEA equations is solved in 
appendix A. 4. Finally in appendix A. 5 the effect of strong Hubbard interaction on the 
A-model and the extended E-model is shown to give respectively the ferromagnetic and 
antiferromagnetic t-J models. 



2 Lattice models based on Shastry's R-matrix 
2.1 The Hubbard model 

The one- dimensional Hubbard model describes spin- 1/2 electrons propagating along a 
lattice with the following nearest-neighbour Hamiltonian 

L L 

H = E = E + 4u , (2.1) 

i=i i=i 

where L is the length of the chain and the kinetic (hopping) and potential terms are 
given by 

V& = ^ (n,,t - J) (n,,, - + J (n,,, - ^) (n,,, - - i . (2.3) 



6 



Here the canonically anticommuting fermionic operators ^ create and Cj cr annihilate 
electrons of spin a =t or a =J, at the j-th site of the lattice. The operator Wj^^ = cj^^Cj^o- 
is the local particle number operator for electrons of spin a at site j, and in what follows 
we also use = rij^^ + n^^^. The coupling constant u is u = U/At, and the hopping 
parameter t is set to 1. We also assume that u is positive. 



The Hamiltonian (2.1 ) is very symmetric. Beyond possessing parity and translational 
invariance it is invariant under a spin 5u(2) and a charge su(2), the generators of which 
are given respectively by 



S+ Y^-^ C+ C— (i~ Q2; 

where 

"it = 4;S\t' ^7 = ^J.tCM> = i(ni-l)- 

Often one considers the more general Hamiltonian 



(2.4) 



(2.5) 



H-/iN-25S" (2.6) 

with N = Ylj=i ~ + L. The chemical potential /i and magnetic field B break 
the charge and spin 5u(2) respectively, leaving only the two u(l) subalgebras unbroken 
as [H, N] = and [H, S^] = 0. Note that 

V"-Evf,^, = X:n„t%,-^N. (2.7) 

i=i i=i 

The Hubbard Hamiltonian is most commonly presented in a form that hides the charge 



su(2) symmetry with /i = — 2u and i? = in (2.6). 

There have been various attempts to extend the Hubbard model by introducing 
correlated hopping and adding various extra terms to the interaction of the Hamiltonian 
(also by introducing next-nearest-neighbour hopping but this shall not be considered 
here). In general one would expect these to explicitly break the symmetries of the 
model. It is also possible however that certain extensions realise different symmetries, 
e.g. [10]- [15], and we see this below in the A- and B- models. 



2.2 Construction 

The Hubbard-Shastry Hamiltonian can be obtained by using Shastry's R-matrix [5] and 
the standard techniques of integrable spin chains [161 S] ■ It was found recently [17] that 
Shastry's R-matrix is equivalent to the 5u(2|2)c-invariant R-matrix [IS] which appears 
as a building block of the AdS^ x scattering S-matrix playing an important role in the 
AdS/CFT spectral problem, see [191 120] fo^' ^ review. Moreover, the su(2|2)c invariance 
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together with the physical unitarity condition is sufficient to fix an R-matrix satisfying 
the Yang-Baxter equation up to a phase [18] . 

The Hilbert space for each site of the one dimensional Hubbard model is a four 
dimensional graded vector space because a site can be either empty or doubly occupied, 
or can have a single electron of either spin. The fundamental representation of su(2|2)c 
is also a four dimensional graded vector space Y , see [Hj for detail. The Hilbert space 
for the model on a lattice of length L is thus naturally identified with 

|0)^|e,,i), 4^4J0) ^ |e,,2), 4^10) ^ |e,,3), 4J0) O |e,-,4) . (2.8) 

The identification is not unique as one can perform canonical transformations on Cj o- , ^. 

A representation of su(2|2)c depends on the values of the central elements. These 
can be conveniently parametrized by variables satisfying 

+ 1-1 , . 

+ — - X = 4zu. (2.9) 

The constraint defines a torus with u a free parameter that characterises its elliptic 
modulus and that we will interpret as a coupling constant. Indeed in the setting where 
one obtains the Hubbard model this is precisely the coupling constant that appears in the 



Hubbard Hamiltonian (2.1 ). The su(2|2)c-invariant R-matrix R(x]'', x^) depends on two 
spectral parameters which from the point of view of the su(2|2)c algebra parametrize the 
two 4-dim fundamental representations upon which the R-matrix acts. We restrict our 
attention to homogeneous spin chains, that the variables of each of the representations 
of su(2|2)c comprising the chain are the same at each site. Then x^ play the role of 
extra coupling constants of the model and the density of the nearest neighbour Hubbard- 
Shastry Hamiltonian is obtained in the standard way through the formula 

Hij+i ~ « Pjj+i9i% j+i(x^, a;^)L±=^± , (2.10) 

where Pjj+i is the permutation matrix, see appendix A.l for precise definitions. 



Any point x^ on the torus can be used to obtain a homogeneous integrable spin chain 
model. In general however the corresponding Hamiltonian will be neither hermitian 
nor parity invariant. One can show by using the properties of the R-matrix that it is 
hermitian if either x^ /x~ or x'^x^ is a phase. In the next section we examine the Bethe 
equations of the models and find that there are just four distinct points on the torus, up 
to symmetries, that give rise to hermitian parity invariant models]^ The Hubbard model 
is obtained when = l/x^. Two of the models have not been singled out before, one 
at x~ = —l/x~^ which we call the A-model, and one at x^ = —x^ which we call the 
B-model|^ Completing the picture is the point x~ = x'^. Here the central extension 

^The same conclusion can be reached by using the transfer matrix and the crossing relation for the 
su(2|2)c-invariant R-matrix. 

•^The A- and B-model Hamiltonians are particular cases of the general Hubbard-Shastry Hamilto- 
nian derived in [6l[7], see eq.(12.229) of [4 . In the parametrization (12.109) used in [4j the A-model 
Hamiltonian (up to a unitary transformation) corresponds to / = ?7r/4, and the B-model corresponds 
to /i = 7r/4. The line x^x^ a phase corresponds to the line Re{fi) — 0, and the line x^/x^ a phase 
corresponds to the line i?e(/i) = 7r/4. 
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of the algebra vanishes and one obtains the su(2|2) spin chain. It is singular from our 
perspective as the coupling u drops out. 

Given a Hamiltonian one can study either the model determined by H or — H. For 
example instead of the Hubbard model one could consider 

L 



H 



attr.H 



E (-T..-M-4uV£^0 . 



(2.11) 



This is the attractive Hubbard model and its physical properties are very different from 
the repulsive case, see e.g. For the A- and B-models, which were constructed not 
on physical grounds but to have rich symmetry and act on the Hubbard lattice, one 
has a choice for the sign of the Hamiltonian. With the hope that they capture some 
interesting physics the sign is chosen so that there is a cost for a doubly occupied site 
relative to a singly occupied site, mirroring Coulomb repulsion. Note however that 
although properties depend strongly on the choice of sign the two cases are closely 



related. This is discussed in appendix A. 3 at the level of the TEA equations. 



Let us comment on the sign of the kinetic term in (2.11). Under the unitary trans- 
formation generated by 

L 

Ui (a) = exp {iaY,3 (%,t + n,u)) (2-12) 



with a an arbitrary real parameter, the kinetic term transforms as {k = j + 1) 



■«a„t 
■k 



(2.13) 



while the interaction is invariant. Thus the Hamiltonian (2.11) is unitary equivalent 
to 



E (T..-M-4uVf,.^,) 



(2.14) 



under Ui(7r). In momentum space X]j=i Tjj+i is diagonahsed as —2 cosp (np_-i-+np j^) 
and the transformation Ui(a) corresponds to a shift p ^ p' = p + a. This is why the 
attractive Hubbard model is conventionally written as (2.14), as states with p' = then 



minimize the kinetic part of the Hamiltonian. In general for a Hamiltonian H we define 
the Hamiltonian H which has the opposite spectrum 

H = -Ul(7r)HUi(7r). (2.15) 

For this reason in what follows we refer to H as the opposite Hamiltonian. 

The direct calculation from the R-matrix at a;~ = 1/a;^ does not give the Hubbard 



Hamiltonian (2.1) but rather one must perform the unitary transformation generated 
by Ui(— 7r/2) to obtain it, which also twists the charge su(2) generators. To obtain the 
Hamiltonians for the A- and B-models (2.16, 2.19) one must apply Ui(7r) and Ui(— 7r/2) 



respectively to the corresponding results from the R-matrix calculation. These twists 



can be understood also at the level of the Bethe equations, see the relations (3.5). 
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2.3 The A- and B-models 



Both the A- and B-models have correlated hopping kinetic terms and we will discuss 
these soon in turn. The interactions of the models are combinations of several common 
terms. In addition to the Hubbard interaction V^, these are 

= = i(n,-l)(n.-l)-i , 

vff = ^(s;s,- + STS+) + s^^s^ 

1 1 



2^ J.I •''-^ ^.-l- -^'1 fs.! '"'-^^ 4 



The charge-charge interaction V'"'^ has the interpretation of nearest neighbour Coulomb 
repulsion. Physically one should think of a site as having charge and so a site 
with one electron on it has neutral charge, a site with electron density less than one is 
positively charged and a site with electron density greater than one is negatively charged. 
The spin-spin interaction V'^'^ is the familiar spin-exchange term of the Heisenberg XXX 
spin chain. Finally the pair hopping interaction V^^ relates to the simultaneous hopping 
of two electrons from one site to a neighbouring site. 

Now we discuss the A- and B-models in turn. The presentation of the Hamiltonians 
is much cleaner if the couphng constant u is re-expressed through 

u = sinh u . 

A-model 

The model at x^x^ = — 1 has the Hamiltonian 

L 



T^fe = + (2.16) 

It is easy to check that it is hermitian but its parity invariance is more subtle. Note that 
the transformation 

CL-j+l,a . ^ CL-j+l,a , (2.17) 

is not a symmetry but is instead equivalent to a replacement u — )■ —u in the Hamiltonian. 



This however is just a unitary transformation. More generally, the Hamiltonian (2.16) 
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Figure 2: The four cases of an electron hopping between sites j and k: i) no electrons of the 
opposite spin on either site, ii) an electron of the opposite spin on site j only, iii) an electron 
of the opposite spin on site k only, iv) electrons of the opposite spin on both sites. 



is unitary equivalent to the following one 

(a) = U| (a) U2 (a) , U2 (a) = exp [ z a ^ n^, -fn^- ; 

r 

2 cosh 2 1/ 



i=i 
2 



Tj^fc(") = - X] [(^^Cfc.a + 4,<7Ci,<T) (1 - n^-,^ - Ilk -a - sin a tanh z/ (rij- - n^.,^)) 



-icosatanhz/(c],^Cfc,^ - c|.^^Cj-^) (n^-^ - iik-a)' 

where a is an arbitrary real parameter. If one chooses a = it one gets 

H^fvr) = with ^ -z/ , 



(2.18) 



and so the parity transformation (2.17) should be supplemented with the change of 



basis generated by U2(vr). Note also that with the choice a = ±tc/2 one removes the 
imaginary term from the Hamiltonian. 

The kinetic part of the Hamiltonian is of a complicated correlated hopping type. 
It is instructive to see how it acts in various cases. For an electron of spin a to hop 
between two sites j and k they must not contain other electrons of that spin but they can 
contain electrons of the opposite spin. We consider individually the four possibilities as 
illustrated in Figure |2] for a mobile spin-up electron. The kinetic density T^f,{a) acting 
on each of the cases takes the form 



i) 



m 



IV 



— ie~*° tanh u cj^^c^^o- + ie*" tanh u 
— ie*" tanh u cj^^Cfc^g. + ie~*" tanh u c[ g., 
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Note the different sign between cases i) and iv), while in the case of free electrons all 
the cases would be the same as case i). This reflects the strong electron correlation in 
the A-model which cannot be reduced to free electrons for any coupling. 

Like all the models we consider, the A-model inherits spin and charge su(2) symme- 
tries from the su(2|2)c structure. The spin su(2) generators take the same form as those 



for the Hubbard model (2.4) whereas the charge 5u(2) generators take a untwisted form. 

The weak coupling z/ — )■ limit of the A-model is not described by free electrons 
but rather the Hamiltonian reduces to the su(2|2) spin chain. This has received much 
study due to an interesting generalisation found in [lOj. A key observation is that 
commutes with H^| q and so one can extend the model, as will be discussed in section 



2.4 In the strong coupling u oo limit only the term 4uV^ in the Hamiltonian survives 



and it thus shares the strict strong coupling limit of the Hubbard model. In section 5^ 
we show that at half-filling and strong coupling the model behaves as a ferromagnetic 
XXX spin chain. The strong coupling limit of the less than half-filled model is analysed 



in appendix A. 5 and is found to be a ferromagnetic t-J model. 



B-model 

The model at /x^ = —1 has the Hamiltonian 



H" = E(Tf;.« + 2'-l"'(v£«-V,%+Vf,,+Vj;«,)). (2.19) 

i=i 

^ - J2 + (l - (1 - (n,,-. - n,,_,)') . 



Let us begin the analysis by examining the correlated hopping kinetic term for the four 
cases illustrated in Figure |2] The model is clearly parity invariant and so cases ii) and 
iii) are identical. The kinetic density T^^(a) acting on each of the cases takes the form 



1), IV) -C'j^^Ck,a - 4,aCj,<7, 



ii) ' iii) - dh^ (S>^^.- + 

In cases i) and iv) hopping is of the standard form, while in cases ii) and iii) its contri- 
bution is proportional to sechu and decreases to zero in the strong coupling limit. Both 
cases i) and iv) preserve the number of doubly occupied sites, while for hopping in cases 
ii) and iii) the number will be either increased or decreased by 1. 

The B-model describes free electrons in the weak coupling z/ — )■ limit. In the strong 
coupling z/ — > oo limit note that [H^ , V^] = 0, as for the weak coupling limit of the 
A-model. This allows one to decouple the coefficient of the Hubbard interaction from 



the other terms while preserving the integrability, see section 2.4 



The B-model possesses a rich symmetry structure. Its spin and charge su(2) gener- 



ators take the same form as for the Hubbard model (2.4). Moreover the model inherits 



a hidden su(2|2) fermionic symmetry for arbitrary coupling constant u from Shastry's 
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Qo,t 


Qo,4. 






Qo,t 








hQ 


2 


2 


-2 


-2 


-2 


-2 


2 


2 




-1 


-1 


-1 


-1 


1 


1 


1 


1 




-1/2 


1/2 


-1/2 


1/2 


1/2 


-1/2 


1/2 


-1/2 



Table 1: Commutators of supercharges of the B-model with H, N and S^. 

R-matrix. There are eight supercharges, four of which Qo^., Qo.a create/annihilate an 
electron with momentum p = 0, and the other four Q]^^^, Qn,o- create/annihilate an 
electron with momentum p = it. Explicitly they are 

L 

Qaa = J^e^°^'(7aC,-,-/3n,-_,c,-,), Q],,, = (Q«,a)^ a = 0,7r, (2.20) 



1 / e*°\ / 1/ 

with 7^ = - ( /3 + —1 and /3 = W tanh - , 

and they satisfy the su(2|2) algebra with central charge cothz/: 

{Qo,;,Qi,t} = -{Qn,i,CiU} = S+ , {Qo,t,Qj,J = -{Q.,t,Qt,J = S-, 

{Qo,t,Qi,t} = S'-^^ + cothz., {Qo,4,QJJ = -S"-77^ + cothz/, 

{Q^,t,Ql,t} = -S" + T7" + cothi^, {Q^,;,q|^ J = S" + T7" + cothi^, 

and all other anti-commutators vanish. The symmetry is hidden in that these generators 
do not commute with the Hamiltonian but rather 

[H^,Q] = /iqQ, [N,Q] = nQQ, [S^ Q] = Q , 

with hq, uq and sq given in Table [ij Thus the symmetry does not lead to an extra 
degeneracy in the spectrum of the model, but it does imply exact relationships between 
the energies of eigenstates. For example if {ip) is an eigenstate with energy then 
the energy of the state Q\ip) is given through H^QIV') = (S^ + /iq)Q|V'). In particular it 
follows that the ground state of is in the kernel of qI,^, Qqi, Qn,t and Qtt,i- 

Let us conclude with two remarks. First, the above fermionic symmetry can be 
generalised naturally to any bipartite lattice in any dimension. This raises the prospect 
of being able to extract some exact information for such lattices. Finally, although the 
symmetry is hidden for H^, it can be realised in part for the Hamiltonian H^— /i N— 2 B 
for appropriately tuned fi and B. To this end it is useful to consider not just the above 
supercharges Q but also products of them with S+, S^, t/+ and rj . 



2.4 Extended models 

The Hubbard interaction term preserves both spin and charge su(2) symmetries. It is 
physically reasonable to add this term with an arbitrary coupling q to the Hamiltonians 
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of the A- and B-models and the opposite models because it allows strong onsite Coulomb 
repulsion of electrons to be introduced. In particular the g — )■ oo limit is studied in 



appendix A. 5 where it is shown that the B-model reduces to the t-J model with the 
coupling interpolating between J = and J = 2 as u changes from to oo. 

In general the extension destroys the integrability of the models but there are special 
values of u for which it does not. In this section we present integrable extensions of the 
weak coupling limit of the A-model and the strong coupling limit of the B-model. We 
borrow an idea from |10] and the model below is indeed the one they studied. The 
key observation is that if for a model [H , V^] = then H and can be diagonalised 
simultaneously and it is natural to add a term 4 q to the Hamiltonian, where g is a 
new coupling constant. This allows one to decouple the term capturing onsite Coulomb 
repulsion from the other terms in the Hamiltonian. Since the number of doubly occupied 
sites is conserved if [H, N] = 0, then q can be thought of as a chemical potential for 
doubly occupied sites. Furthermore if the original model is integrable and commutes 
with the transfer matrix then the resulting model is integrable too. This is the case for 
the limits of the A- and B-models that we consider here. 

The extended Hamiltonians for the weak coupling limit of the A-model and its op- 



posite model (see equation (2.15)) are 

H^°(fl) = E(Tj;°,, + 2((20 + l)V£,, + V5:f,,-Vf,,+Vj;f,O),(2.21) 

H^°(fl) = X:(Tj:?.i + 2((20-l)V£,,-V$^,+Vj^,-Vj-^O)>(2.22) 



The Hamiltonian H'^°(s) is the one for the EKS model [TO] (with q = U/A), and eq.( |2.22[ ) 
shows that for g < 1/2 it should be thought of as a model with attractive on-site Coulomb 
interaction. The extended Hamiltonians for the strong coupling limit of the B-model 
and its opposite model are 

H"-(0) = E(TL>i + 2((20 + l)Vf^.^,-V$;f^,+V^;f^, + Vj:f^O),(2.23) 
H^-(0) = E(T,^.-i + 2((2s-l)Vf^.^,+V$;f^,-V^?f^,-Vj;f^0),^^^^^^ 



= - E (S- -C'^.- + (1 - (n„ - ^k,-.?) . 

(7 

Similar to the B-model, both extended models have hidden fermionic symmetry. 
Since the extended A-model with g = is the su(2|2) spin chain its Hamiltonian com- 
mutes with the supercharges (TU] which create and annihilate electrons on singly occupied 
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and doubly occupied sites 

L L 

Qs,o- = ^ ] (^jiO' ~ ' Qrf.o- = ^ ] Hj,-o-Cj,o- , (2.25) 

and their hermitian conjugates. Noting that 

[v^ , Q.,.] = ^Q.,. , [v^ , Qd,a] = -^Qd,. , (2.26) 

one finds 

[H^"(0) , Q.,<x] = 20 Q.,. , [H^'Xfl) , Qd,a] = -20 Qd,a ■ (2.27) 
Taking the hmit u — )■ oo for the supercharges ( |2.20[ ) of the B-model one finds similarly 

L L 

Qo,a = X^(Cj> - Hj^^Cj^^) , Q^_^ = ^ (-l)^nj _^Cj- <^ , (2.28) 
j=i j=i 

and 

[H^°°(0),Qo,<x] =2(0 + 1) Qo,., [H^°°(0),Q.,.] = -2(0 + 1) Q,,,. (2.29) 

The Hamiltonians and symmetries of these models obviously look different but, as 
we will see in section [7| for even L they have the same Bethe equations. Analyzing 
their dispersion relations one finds that the Hamiltonians H^°°(0) and H^o(0 + 1) have 
the same spectrum and therefore they should be unitary equivalent. The corresponding 
unitary transformation is given by 

U = exp TT ^ j rij, ^n^,; j , (2.30) 

and it shifts the momenta of doubly occupied sites by vr while leaving the momenta of 
singly occupied sites unchanged. Correspondingly it untwists the supercharges Q^^.o- of 
the extended B-model but leave unchanged Qo,cr- 

3 Bethe equations 



In the previous section a family of models with Hamiltonian (2.10) were introduced. 
These were constructed from the su(2|2)c invariant Shastry R-matrix using the quantum 
inverse scattering method. The models have been diagonalised for arbitrary using 
the algebraic Bethe ansatz [2T| [22l [20] . For the purpose of identifying the momentum 
and dispersion relation, the relevant part of the eigenvalue of the transfer matrix t for 
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a state on a periodic one-dimensional chain of length L with N charge excitations over 
the pseudo-vacuum (|0))^ is given by, see e.g. [171 [23] 




t=\\e^'^^-^^-^J^, (3.1) 



k=i ^'^ 



where G M is an arbitrary twist which should satisfy the condition exp(i0L) = 1 for a 
periodic spin chain. The yk are the Bethe roots for charge excitations. The momentum 
is then defined through 




em- = e«0il^ J — ^ -vr<pfc<7r, (3.2) 



as the transfer matrix is the shift operator and so its eigenvalue is naturally identified 
with e*^'--^''. The roots yk are determined from the Bethe equations 

M , . 



,^P.L^TJ -k-u.,^.u ^ fc=l,...,iV <L, 

Vk — Wj — iu 

(3.3) 



^v^k^v^^ ^ Wk-w. + 2^u ^ fc = l,...,M<iV/2, 
Wk — Vi — iu Wk — Wj — 2iu 

where w are the Bethe roots for spin excitations]^ and the v are related to the charge 
Bethe roots via v = ^{y + 1 / y) . 

As stated in the previous section the models will only be hermitian if either /x^ 
or x~^x~ is a phase. The requirement of parity invariance is more stringent]^ At the 
level of the Bethe equations it amounts to a symmetry under p — )■ —p. There are two 
distinct ways this can be realised 

p — )■ —p, y — !■ —y, w — )■ —w , 

1 (3.4) 
p ^ —p, y — )■ — , w ^ —w . 

y 



However p and y are related through (3.2) and so together with equation (2.9) these 
transformations strongly constrain x^ and (p. One finds that there are four possible 
solutions leading to four distinct parity invariant models, each of which satisfy the 
requirement to be hermitia]:]0 

(H) x+ = l/x- = oo, 0=2z^7r, 

(A) x+ = -1/x- = i{u + Vl + u2), (f) = n7c, (3.5) 

(B) x+ = -X- = i{u+Vl + u^), 0=^7r, n = 0,l. 



■*Our notations are related to the ones adopted in [3] as pj -n- kj, wj O Aj. 
^We thank Niklas Beisert for a valuable discussion of this issue. 

^The fourth solution = x~ = oo corresponds to the su(2|2) spin chain but this point is highly 



degenerate and to get the correct Bethe equations for the model from (3.3) one should send the roots 
Vk and Wk to infinity in a proper way. Technically it is similar to the way we get the Bethe equations 
for the Bu(2|2) spin chain from the weak coupling limit of the A-model in sectionj?] 
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In what follows we choose (/)h = — f , 0a = tt and 0b = — | wh ich lead to the minus sign 
in front of the cosp term in the dispersion relations, see (3.6). 

Next we wish to obtain the dispersion relation of these models. To proceed one 
observes how the Hamiltonian acts on one particle states. The result can be obtained 
from the explicit forms of the A- and B-model Hamiltonians (2.16 2.19), but it is 
of interest to obtain it in the present context also. The transfer matrix is related to 
the Hamiltonian through H ~ it~^dt, and so its eigenvalue for a single particle, the 



dispersion relation, follows straightforwardly from the transfer matrix eigenvalue (3.1) 
with = 1. For the four models, using the same normalisation as for the Hamiltonians 
of the previous section, one finds 



(H) 


m = 


— 2cos(p) — 2u. 


(A) 


m = 


-2cos(p) -2^1 + u2 


(B) 


m = 


—2 cos(p). 



(3.6) 



This clearly agrees with what one expects from the explicit forms of the Hamiltonians. 

Let us conclude by remarking how the transformation between a model and its op- 
posite, given by equation (2.15), can be understood in the present context. Shifting 



p — )■ p + TT, which is equivalent to choosing the other value for n in (3.5), and changing 



the overall sign of the dispersion relation (3.6), provides the dispersion relation for the 
opposite model. 



4 Thermodynamic limit 

Let us turn our attention to the physical properties of the models introduced above. We 
will examine the thermodynamic limit, L — )■ cx3, in the grand canonical ensemble. The 
equilibrium thermodynamic properties follow from the minimization of the free energy 



/(/i,i?,T) = -^log(tr 



exp 



H-/iN-25S' 



(4.1) 



Here ix is the chemical potential, B is the magnetic field, and T is the temperature. 

The magnetization m per site, particle density rzc, specific heat capacity c, and the 
spin Xs and charge Xc susceptibilities are calculated in the standard way 

Idf df rpd'^f 9 m dric 

The free energy obeys the following symmetries 

/(/i, -B, T) = fifi, B, T) , /(-/i, S, T) = /(^, B, T) + 2/i . (4.3) 

The first follows from the invariance of the Hamiltonians of the models with respect 
to spin reversal. To see the second note that N — t- 2L — N under the particle-hole 
transformation cj^ ^ (— l)-'Cj _a. for the Hubbard and B-models, and cj^^ ^ ^j-a for 
the A-model, which is also a symmetry of the Hamiltonians. 
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4.1 String hypothesis 

The free energy is calculated in the thermodynamic limit using the Bethe ansatz and 
making a string hypothesis [21] for the behaviour of the roots in this limit. 

To describe the string hypothesis and TEA equations it is convenient to introduce 
the following functions 



xa{v) = v + v^1- ^, \xA{v)\>l,veC, (4.4) 

with a cut = (—1, 1), and 

xb{v) = v + iVl - , lm{xBiv)) > , V e C , (4.5) 

with a cut = (— oo, —1) U (1, oo). Let us also define X^'^ and X^'^ to be respectively 
the upper and lower edges of X^'^. For values of v on the cuts we define x{v) = x{v + iO). 
Both functions solve the equation 

x{v) + = 2v , 

X[V) 

and the parameters x^ of the A- and B-models can be uniformly written as x~^ = x{iu), 
x~ = l/x{—iu) where x = xa and x = xb for the A- and the B-model, respectively!^ 

Since y + 1/y = 2v for any given v one has two possible corresponding y-roots and 
they can be parametrized as 

y+ = x{v) , y_ = ^— , (4.6) 

X[V) 

so that the set of complex y-roots is divided into these two subsets. Here and in what 
follows x{v) = xa{v) for the Hubbard and A-models, and x{v) = xb{v) for the B-model. 
The two distinct types of |/-roots y± obviously lead to two types p±{v) and S±{v) of the 



momentum p and the energy S as functions of v, see (3.2) and (3.6) respectively. 



The string hypothesis, see [211 |25l HI [26] , states that each root of the Bethe equations 
in the thermodynamic limit is a member of the following types of Bethe strings: 

1. y-particle: a single charge spin-up excitation parametrized by its real momentum 
p. The reality of p then implies that 



IvI = 1 , If I < 1 , f G R for Hubbard and A — models , 

' ' ' (4.7) 

y G R, |f I > 1 , fGR for B — model. 



^ The parameters of the two famihcs of Hermitian (but in general not parity-invariant) Hubbard- 
Shastry models can be written as a;+ = x{X + iu), = l/x{X — iu) where A is an arbitrary real 
number. For the models with x^x~ being a phase x — xa, and for the ones with x'^ /x~ being a phase 
X = Xb- Then the Hubbard and the su(2|2) spin chain model's parameters correspond to A = +00. 
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Since v is on the cut of the corresponding function x its y± roots can be also 
written as 

y± = x{v ± iO) . (4.8) 

This formula implies that for single charge excitations p±{v) and S±{v) can be 
thought of as the values of the functions p{v) = p+{v) and S{v) = S+{v) on the 
upper and lower edges of the cut of x{v). One finds for all the models that moving 
around the cut of x{v) in the counter-clockwise direction increases the momentum 
p, and, therefore, dp-/dv > 0, dp+/dv < for v being on the cut. 

2. Ml^uJ-string: a spinless charge 2M bound state composed of 2M roots yj and M 
roots Wj and parametrised by v e R 

Vj = v + {M + 2-2j)ui, v_j = V - {M + 2-2j)ui, (4.9) 
Wj = v + {M + l-2j)ui, j = l,...,M, (4.10) 

where the roots yj are expressed through vj as 

yj = x{vj) , y_j = l/x{v_j) , j = 1, . . . , M . (4.11) 

As a result, the momentum and the dispersion relation for the Ml^w-string are 

M M 

PM\vwiv) = b+ (vj) + P- {v-j)] , Sm\vw{v) = Yl [^+ (^^■) + ^- (^-:')] ■ 

3=1 3=1 

(4.12) 

One can check for all the models that dpM\vw/dv < for any real v. 

3. M|w-string: a spin — M bound state composed of M roots wj and parametrised 
by w e R 

Wj^w + {M + l-2j)ui, j = l,...,M. (4.13) 
This family includes the l|i/;-string which has a single real root w. 

The string hypothesis for the Hubbard and A-models differs from the one for the B-model 
only by the location of a root y of a single charge excitation. This is accommodated by 
the functions xa,b{v), and in terms of these all the formulae look the same. 

4.2 Free energy and TBA equations 

In the thermodynamic limit we introduce the densities pk and pk of particles and 
holes, and the Y-functions = Pk/pk and dressed energies = TlogYfc, where 
k = ± , M\vw , M\w , M > 1. Since, as was discussed above, the roots y± are equal 
to the values of x{v) on the upper and lower edges of its cut, the Kt-functions satisfy 
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the same property and can be viewed as the values of the function Y{v) on the upper 
and lower edges of its square-root branch cut coinciding with the cut of x{v) 

Y+{v) = Y{v + iO) , Y^{v) = Y{v -iO) . (4.14) 

The derivation of the free energy and TBA equations follows a textbook route [1], 



see appendix A. 2 for some detail. The free energy is given by 



1 — - ^ dp T ^ ^ f dp dpi|^^ 



/ = + -(f - f,„„ -)« d§ - ^ (1 + V) « K§ - - ^ ) (4.15) 



T dv 

■^iog(i + yi|.)vc.®^, 



where the symbol * denotes the following "convolution" 

POO 

g'kh= / dt g{u,t)h{t,v) 



(4.16) 



and ® denotes the contour integral in the counter-clockwise direction around the branch 
cut of Y{v) or, equivalently, of £{v) and p{v) located on the real axes. Explicitly for the 
Hubbard and A-models one has 



g ® h 



dt {g_{u,t)h-{t,v) - g+{u,t)h+{t,v)) = g-^h- - g+^h+ 



(4.17) 



while for the B-model 



g®h= I dt {g^{u,t)h_{t,v) - g+{u,t)h+{t,v)) = g_^h_ - g+^h+, (4.18) 
'\t\>i 



where 



g±{u^ t) = g{u, t ±iO) , h±{t, v) = h{t ± iO, v) 



(4.19) 



and ^ and denote convolutions with the integration over |t| < 1 and |t| > 1 re- 
spectively. If g (or h) is a function of a single variable then one just drops u (or v or 
both) in (4.16 4.18), e.g. g = g(t) then g -k h = dt g(t)h(t,v). However, ii g or h 
is a kernel denned through a function of one variable then it should be understood as 



g{t,v) = g{t — v). Finally in (4.15) s denotes a kernel defined in (A. 1.6) 



The F-functions in (4.15) are determined from the simplified TBA equations 

logy(f) = ^ +log ' -kS, 

log Ya//|,;«;(^^) = /mat log (l + YatI^^) *S - 5m1 log(l + F)®S , 
\ogYM\w{v) = /M7vlog(l + Y)v|^„)*S-5Mllog(l+ ^)®S, 



(4.20) 
(4.21) 
(4.22) 
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which are the conditions that the free energy is indeed minimised. Here Imn = 5m+i,n + 
8m-i,n , M >2 and Iin = S2n- Note that eq.(4.20) imphes the following simple relation 
between Y+ functions 



log 



T 



(4.23) 



following from the analyticity of the kernel s{t — v) for real v. 

In addition the TBA equations are supplemented with the asymptotics for large M 



lim 



log Ymk 

M 



2/i 



lim 

M^oo 



log Ya/ K 
M 



2B 



(4.24) 



which follow from the canonical TBA equations, see appendix A. 2 for the details. Note 
that the chemical potential /i and magnetic field B enter the simplified TBA equations 
only via the large M asymptotics of Ym\vw and Ym\w 



It often happens that the sign of the final terms in equations (4.21), (4.22) play an 
important role in understanding the physics of the models. Through equation (4.23) 
these signs follow from the relative magnitude of and Let us thus note that 



-4Vl - v'^ < 



V? + f ^ 



> 



> 



for the Hubbard model, 
for the A model, 

for the B model. 



(4.25) 



In the case of the Hubbard model one can easily check that the TBA equations and 
the free energy match exactly the ones in j^. To be precise one finds the relations 



eo = u + —{£ -£i\i 
Zn 



-k s 



dp 



1 

2^ 



dp 



— Sic 



dp 
dv ' 

dpl\^ 



27T 



dp 
dv 



(4.26) 
(4.27) 



dv 271 V dv dv 
where cq, cq and po are the ground state energy per site, the density of w-strings, and the 
density of roots of single charge excitations for the half filled repulsive Hubbard model, 
respectively, see eqs.(5.69) and (5.70) of [1]. 

Finally it is worth pointing out that, as follows from the TBA equations, the Y- 
functions have infinitely-many square-root branch points on the f-plane. To be precise, 
Y{v) has branch points located at v = ±1 + 2kui , k E Z, while Ym\vw and Ym\w have 
branch points located at v = ±1 + M ui + 2kui and v = ±1 — M ui — 2kui , fcGN. 



5 Various limits 



The free energy (4.15) together with the TBA equations (4.20)-(4.22) and their large M 



asymptotics (4.24) contain all the information about the equilibrium states of the Hub- 



bard, A- and B- models in the thermodynamic limit. In general one needs to numerically 
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solve the coupled integral TBA equations to extract this information. In various limits 
however the equations simplify and one can make analytic progress. 

We first consider the weak and strong coupling limits. As the kernel s{v) has explicit 
dependence on u, it is useful to rescale the variable v to v/u and introduce a notation 
for functions that take a rescaled argument 



yk{v) = Yk{uv) , 5{v) 



4 cosh f ' 



1 



M 



(5.1) 



Here we have also altered s{v) and Km{v) by a factor of u for convenience as they 
generally appear together with dv. Note that after rescaling the variable in the TBA 
equations they take the same form with Y-functions replaced Yk by yu and s by s. 

The limits of strong magnetic field, large chemical potential and infinite temperature 
are also discussed. Let us here introduce the convenient notation 

/3 = L /3c /3. = ^. (5.2) 



T 



T 



5.1 Weak coupling limit u — > 
Hubbard and B model: free fermions 



From the Hamiltonians of the Hubbard model (2.1) and the B model (2.19) it is clear 



that they reduce to free fermions in the weak coupling u — ?■ limit. Here the limit is 
analysed on the level of the TBA equations. We find as expected that TBA equations 
can be solved explicitly and thus a closed expression for the free energy can be obtained. 
Let us consider the u — )■ limit without rescaling the variable f |^ First note that 



limu_>.os(v) = Next for v on the intervals = (—1,1) and = (— oo, 

(1, oo) for the Hubbard and B-models respectively, one finds the following limits 



U-5>0 

hm — — (' 
u-s>0 dv 



A 



lim£^_ 

u^O 

2 

A ' u"-To dv 



-A, 



lim £i 

U-!>0 



lim 



dp- 



2 

A' 



lim 



\vw{V) -- 

dpi|t 



2A, 



-1)U 



(5.3) 



dv 



4 

'a 



where for the Hubbard model A = 2a/1 — v'^ and for the B-model A = —2^1 - 
Then for these values of v the TBA equations become a set of algebraic equations 



iogy+ 

^OgYM\w 



Y 



1, 1 

- log — 

2 ^ 1 + Fii 



log 



2A 



Imn log (l + Yj 



N\vw^ 



+ 2^Milog^^ 



IMN 



lo. 



1 



1 + 



g (l + YN\yj) + -Smi log ^ 

^ J- + \7- 



(5.4) 
(5.5) 

(5.6) 



^It is worth pointing out that if in the Hubbard case one would rescale v then the naive u — > limit 
would lead to a constant solution of the TBA equations and to divergent integrals in the free energy 



expression (4.15), while in the B case one would loose all the information about y±-functions. 
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Note that to compute the free energy (4.15) one needs to know values of Y-f unctions on 
and only. The general solution to these algebraic TBA equations is well-known 



and given by (A.4.53). 



In the limit u — )■ the free energy (4.15) simplifies dramatically 



/ = /dt;-- /^iog((i + y_)v/rTiv) . 

Jx JX ^ 



and by using (A.4.53) one finds 



/ = -Ty ^ log (^(1 + e T )(i + e T ) 
This is the expected result for free electrons [25] . 



A-model: 5u(2|2) spin chain 

The weak coupling limit of the A-model is studied in detail in section [7] in the more 
general context of extended models. Here we point out that the limit can be taken at 
the level of the TBA equations. However if one tries to take the limit keeping v fixed 
then one finds that most of the functions in (5.3) exhibit singular behaviour at w = 0. 
This can be resolved by rescaling v to vju. Note that then the cut (—1,1) — > (—00,00) 
as u — )■ and thus y±{y) can be considered as distinct functions defined on M. In this 



way the TBA equations (7.3) with = can be obtained without going back to the 



Bethe ansatz for the su(2|2) spin chain. 



5.2 Strong coupling limit u — > 00 
Hubbard and A-model: su(2) spin chain 

Here we consider the models in the limit of infinite coupling. The analysis is quite similar 
for the Hubbard and A-models and so these are examined together. 

It is well known that the strong coupling limit of the less than half-filled Hubbard 



model is the t- J model (see appendix A. 5 for a discussion), which at half-filling takes 



the form of the antiferromagnetic XXX spin chain 



H = -2u - J (s,.S,+i - ^) - 25 S% 

i=i ^ ^ 



(5.7) 



where for later convenience we use J = — - . The free energy and TBA equations for this 
Hamiltonian are known to be [25j: 



/ = -2u + Jlog2-T 



dv log(l + 3^1(1;)) 5(t;) 



log3^Af(^^) =/AfJvlog (1 + 3^jv) *S (1;) + 5 



M\- 



27tJ3{v) 



lim 



logy 



M 



2B 



{5.t 
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The Hubbard free energy and TEA equations reduce to these as u — )■ oo. 

Likewise we will see that in the strong coupling limit the free energy and TEA 
equations of the A-model reduce to (5.8) with J = ^. Thus it is natural to conclude 



that the strong coupling limit of the half-filled A-model is the ferromagnetic XXX spin 

1 

u ■ 



chain with Hamiltonian given by (5.7) with J 



First as u oo note that £^ ^ E"^ ^ -2u, £^^^ -k s ^ ^ and £{ 



i\vw 



log 2 



u 



Then from TEA equation (4.21) it follows that F — )■ 0. Furthermore 



log(l + -)®s = log ^ ^ 1 



1 + ^. I,, 



£)is 



y+ 



by equation (4.23). The free energy and TEA equations for the Hubbard and A-models 



thus simplify to 



/ = -/i - 2u + - £ log (1 + Fii^) ^ s i 



dp 
dv ' 



IokF, 



M\w[ 



Imn log (1 + Yni^,) *s-^-^i£+- £^)ks 



(5.9) 
(5.10) 



with K = — 1 for the Hubbard model and k = 1 for the A-model. Since n enters only 
the first term of the free energy it follows that the electron density is one, i.e. it is at 
half-filling. These resemble quite closely the corresponding equations for the XXX spin 



chain (5.8). Eefore the exact equivalence is shown let us remark that the sign of the final 



term in (5.10) gives an indication of whether one should expect an antiferromagnetic 



or ferromagnetic spin chain. From (4.25) we thus expect the Hubbard model to exhibit 



antiferromagnetic behaviour and the A-model to exhibit ferromagnetic behaviour. 
Next note the strong coupling limits 



j\t{£+{t)-£^{tMt,v)^[ 



27r s(v) 



for the Hubbard model, 



-27rs(f) for the A-model. 



Rescaling the kernel s and Y- functions as in (5.1), eq.(5.10) takes the form 

27r5(f) 



yM\wiv) = Imn log (l + 3^Ar|^) -ks{v) + k5_ 



AIl- 



uT 



To take into account the rescaling in the final term of equation (5.9) note 



log (1 + Fii^) ^ a ®g = JZ log (1 + yiUt)) /_\ dvs{t, uv) {^{v) - 



dv 



J^^ dvs{t, uv) 



dp- 
dv 



dv 



large u 



> 271 s{t) . 



Hence under the identification yM\w = 3^a/ the free energy and TEA equations for the 
Hubbard and A models do indeed reduce to those for the XXX spin chain (5.8) with 

J = -K 



J = ^ respectively. 
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Let us comment on the attractive versions of the models. 

u ' l\vw u 



Here £^ ^ 2u, 



^i\vw * s ^ ^i\vw ^ s ^ ^ and so the TEA equation ( |4.2l| ) imply Y oo. It 

follows that the free energy (using the form (A. 2. 44)) and TEA equations simplify to 



/ 



logYAfK 



= Imn log (1 + YN\,y,) ^s+^-^{£+-£-)is 



with K = —1 for the attractive Hubbard model and k = 1 for the attractive A-model. 
Since there is no B dependence in the equations the magnetisation vanishes. As the 



inequalities in (4.25) are reversed in going to the opposite models one sees that in the 



strong coupling limit the attractive Hubbard model is an antiferromagnetic XXX chain 
while the attractive A-model is a ferromagnetic XXX chain. Note that these are charge 
su(2) chains, governed by 77, as opposed to spin chains. For the strongly attractive A- 
model one thus expects spontaneous symmetry breaking at T = with the model being 
empty for /i < and having all sites doubly occupied for /i > 0. 



B-model 

The strong coupling limit of the E-model is studied in detail in section [7j The analysis 
is similar to that for the weak coupling limit of the A-model above. Note that if v is 
rescaled to v /\x then the cut (— cx3, —1) U (1, 00) — )■ (— cxo, 00) as u — )■ 00. Again y±{v) 
can be considered as distinct functions defined on M, and the TEA equations ( |7.3 ) for 
the extended E-model with = can be found. 



5.3 Strong magnetic field limit 5 — > 00 

In the limit 5 — t- 00 the Y-f unctions for M\w strings diverge. This can be seen from the 



canonical TEA equations (A.2.24)-(A.2.26). Note that the final term on the right hand 



side of equation (A. 2. 26) is finite as -B — )■ 00 because y — )■ from (A. 2. 24) and thus 



log(l + -)®Km -£®Km , 



using (4.23). Therefore 



and the free energy simplifies to 

f = -li-B 



27r 



{£ — £i\i 



-kS] 



dp 
dv 



The electron density is one and the magnetisation per site is one-half, and thus the 
ground state is full of spin-up electrons as expected. 
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5.4 Large chemical potential limit /i — oo 

The analysis of the hmit /i — t- — oo is similar to the limit B ^ oo. Here the Y- 



functions for M|fw-strings diverge as follows from the canonical TEA equation (A. 2. 25) 
and become 

J- M\vw ^ ^ 1 



and from (A. 2. 44) one finds the free energy 

1 



/ 



27r 



-k S I 



dp 

dv 



Thus the ground state is empty of all charge excitations. This is as expected as a large 
chemical potential means there is a large energy cost for having a particle in the system. 



5.5 Infinite temperature limit 

In the infinite temperature limit one expects many of the details of a system to become 
irrelevant. In particular we find that the Hubbard, A- and B-models become identical. 
In the limit T — )■ oo, with Pc and Ps fixed, the driving term of the TEA equations 
(4.20)-(4.22) becomes negligible. The Y-functions are thus constants and = Y_ = Y. 



The TEA equations become recursion relations with the solution [1] 



Y 



cosh (3c 
cosh f3s 



Yf 



M\vw 



sinh(M + l)l3c 



sinh (3c 



1 , Ym\i 



sinh(M + 
sinh (3s 



and the free energy becomes 



/ = -/i - riog(i + Y)-^ iog(i + 



-Tlog {2e'^^ cosh/3, + e^^^ + l) . 



In the strict infinite temperature limit /3c and (3s — )■ 0. Thus / — )■ —Tlog 4 as 
expected for a model with four degrees of freedom per site. 



6 Zero temperature limit 

The T — )■ limit is of great interest as the absence of thermal fluctuations allows one to 
obtain a clear picture of the phase diagram. The TEA equations simplify dramatically 
and it is possible to study them analytically in certain regimes. Indeed in this limit it 
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is useful to work with the canonical TBA equations (A.2.24)-(A.2.26| 



\OgYy 



log Ymk 



£y — jJ, — B 



T 

£m\vw ■ 



log 



1 + 



1 + 



-2Mu , / 1 



Km + log 1 + 



Yr 



N\vw 



i<K 



NM, 



logy, 



2MB 



M\w 



log 1 + 



1 



Km + log 1 + 



Y^ 



i.K 



NM- 



y / \ ^ N\w , 

These equations are the direct result of the minimisation of the free energy and it is from 



them that the simplified TBA equations are derived, see appendix A. 2 They contain /x 



and B explicitly rather than through the asymptotics (4.24) which will not be relevant 



for the discussion of zero temperature. Furthermore, the Y-functions will vanish and 
diverge as T — > and so we will work with the dressed energies which remain finite 



efc = TlogYfc, 
where k labels the string type. Note then that 



efc < 
Cfc > 



limT^o Yk 
limr^o Yi, 



0, 



(6.1) 



(6.2) 



Recalling = pk/ pki the ratio of the density of holes to the density of particles, it 
follows that ek{v) < means that there are no holes of type k with spectral parameter 
V in the ground-state whereas ekiv) > means that there are no particles of type k 
with spectral parameter v in the ground-state. From the canonical TBA equations it is 
seen that a Y-function does not contribute if Yk{v) > 1 for all v, i.e. that there are no 
particles of that type. 

It is possible to rule out a lot of string types by momentarily considering the simplified 



TBA equations (4.20)-(4.22). Note that 



Y, 



M\vw 



> 1 and Ym\w > 1 for M > 2 , 



as the right hand sides of the equations for these strings are strictly positive. The sign 
of the extra term in the TBA equations for Yi\y^ and Yii^^ follows from the inequalities 
(4.25). Thus Yi\yy^ > 1 for the Hubbard model as is well known, whereas for the A- and 
B-models one has that Fii^ > 1. 

With this insight the canonical TBA equations simplify dramatically at zero tem- 
perature. From now on we consider only the cases of the A- and B-models. Writing the 
equations in terms of the dressed energies they take the form 



£ . 



fi-B 
, - 2/i4 



ei|t 



K2 



(6.3) 



^In this section we use the subscript y to label quantities related to y-particles. 
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where Kn = K2 was used. Here a subscript notational convention has been adopted 
for the convolutions to signify that the integrals should only be taken over the region 



for which < 0, as for > one has that — )■ 00 as in (6.2) and the contribution 
vanishes. Here Qk is the interval defined by 

ekiy) < for V e Qk- 

We see from the definition of Y-functions that Qi\vw ^ and Qy = U Q^, where 
Q± C X^'^ (recall that I^'^ are the upper and lower edges of the cuts of x^'^{v)). 

The phase diagram is a description of the constituent particles of the ground state as 
a function of fi and B. Recall that the Bethe ansatz deals with states in the regime /i < 
and B > and so it is for this quadrant that the TEA equations provide information. 
The phase diagram for the other three quadrants follows straightforwardly from the 



symmetries (4.3) of the free energy. 

The ground state can consist of only y-particles (spin-up electrons) and l\vw strings 
(bound states of a spin-up electron and a spin-down electron). There are five phases 
that will be relevant to our discussion: 

I. Empty band: no y-particles, no l|f w-strings. 

II. Partially filled and spin polarised band: some y-particles, no l|f w-strings. 

HI. Half filled and spin polarised band: saturated with y-particles, no l|fw-strings. 

IV. Partially filled and partially spin polarised band: some ^/-particles, some l\vw- 
strings. 

V. Half filled and partially spin polarised band: some ^/-particles, saturated with 
l|fw-strings. 

Further phases such as partially filled and spin neutral are not identified as such will 
correspond to points and lines rather than regions of the phase diagram. 

Boundaries between phases are determined as follows: the ground state begins to 
contain a string of type k when and B are tuned so that mm^ek{v) = 0, and the 
ground state becomes saturated with a string of type k when fi and B are tuned so that 
maxt, efc(f ) = 0. 

A-model 

In the A-model there are no l|fw-strings in the zero temperature ground state. This 



can be seen from a numerical study of the TBA equations (6.3). One finds that ei|„^ < 
is not consistent with B > and so there can be no l|fw-strings. Thus the model can 
only exhibit phases I, II and III, and the relevant TBA equation is 

ey = £y-f,-B. (6.4) 
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Recalling the dispersion relation from (3.6), we see that the ground state is empty for 



11 < —2 — 2 a/1 + — B, then partially filled up to the line n = 2 — 2\/l + — B, and 
saturated with y-particles thereafter. The phase diagram for the A-model for /i < 0, 
5 > is presented in Figure [1} Note that i? = is a singular point at T = arising 
from spontaneous breaking of su(2)s. 

B-model 

In the B-model there may exist both y-particles and l|fw-strings in the ground state. 
Let us determine the phase diagram on a case by case basis. 

First consider the question of when l|fw-strings begin to appear in the ground state, 
this determines the boundary between phase IV and phases I, II and III. When there are 
no l|f i(7-strings then ei|„^(f ) > for all v. Thus the convolutions ei|^,^ do not contribute 



to the right hand side of the TBA equations (6.3), and they become 



ey = Sy- fi- B, ei|^^ = - 2/i + Ki. (6.5) 

The phase boundary is determined by the condition that min„ei|^^(f) = 0. Note that 
ey(v) is even and a monotonically increasing function of v around the cut (l,oo) from 
the minimum of Sy{v) at f = oo + to the maximum of Sy{v) at v = oo — iO, that is 
e+{v) and e_(w) are decreasing and increasing for positive v respectively. Also ei|„^(f) 
is even and a monotonically decreasing function for positive v. Thus the minimum 
of ei\vw{v) is at f = ±oo. By evaluating the dressed energies at ±oo one can see that 
they have finite asymptotics, denoted by e^. To do so requires the knowledge that for 
a function f{v) which asymptotes to /°° one has 



/•oo /"OO 

lim / dt f{t) Ki{t,v) = lim / dt f{t + v) Ki{t,0) = f°° , 



for any finite A. These can be used to determine the phase boundaries for the following 
cases i) Qy = 0, the boundary between phases I and IV, ii) C U X^, a proper 
subset, the boundary between phases II and IV, iii) Qy = X^UX^, the boundary between 
phases III and IV: 



i) Evaluated at f = oo the TBA equations ( |6.5[ ) reduce to the linear equations 

X) 

L|dui 



e^ = -2-/i-5, e^ = 2-f,-B, e??,^ = -4 - 2/i. 



From the conditions e'^ > and e^^^ = it follows that the boundary of phases 
I and IV is the point (/x, B) = (-2, 0). 

^°This second property is seen via the identity Siivw — '~£y®Ki — —ey®Ki, where the first equahty 
follows from shifting the contours of integration and the second from the fact that the constant terms in 
Cy cancel when integrated over both edges of the cut. This implies that = — ej,®Q^ Ki, where 
Qy is the complement of Qy in the contour Ul5 around the cut, and it follows that the derivative 
of ^iivwiv) is negative for v £ (0,oo). 
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ii) Evaluated at f = oo the TBA equations (6.5) reduce to the hnear equations 



e^ = -2-/i-i?, e^ = 2-fi-B, ejj^^^ = -4 - 2/i - e^. 

From the condition e^^^ = it follows that the boundary of phases II and IV is 
the line B = 2 + ^. 



iii) Evaluated at f = oo the TBA equations (6.5) reduce to the linear equations 

e^ = -2-ix-B, e-=2-/i-i?, e?f,^ = -4 - 2/i + 4. 

From the conditions e°° < and e^^^ = it follows that the boundary of phases 
III and IV is the point (/i, B) = (0, 2). 

Next let us determine the phase boundaries for y-particles. Neglecting for a moment 
the contribution from l|fw-strings, the relevant TBA equation takes the form 



Then recalling the dispersion relation (3.6), in a similar fashion to the A-model, one can 
see that the line /z = —2 — B would separate phases I and II and the line ^ = 2 — B would 
separate phases II and III. Note that neither of these lines enter the region, located in the 
previous paragraph, for which l|fw-strings appear in the ground state. Thus neglecting 
the l|f t/7-strings in the determination of the y-particle phase boundaries turns out to be 
self-consistent. 

Finally consider when the ground state becomes half-filled. Let us begin by es- 
tablishing that this happens when Qi\yw = The number of y-particles per site is 
ny = I p+ +1 VcQ_ p_, the number of l|t;w-strings per site is ni|^^ = 1 *Qi|„„ pi\vw, 
and then the total number of electrons per site is ric = Uy + 2ni\^^. Consider the zero 
temperature equation for the density of l|fw-strings 

1 dpi|^^ 

Note that Qi\vw =1^ implies 

ni\vw = ^ Pl\vw = l-l-kKiicQ^p^ -l'kKiicQ_p_ -licK2'kpi\y^ 

1 Tiy '^1\VW^ 

and hence ric = I and the ground state is half-filled. It is clear that Qi\vw C M, a proper 
subset, implies ric < 1. Thus the boundary between phases IV and V is determined by 
setting Qi\vw = and varying Qy. First consider the case Qy = 0. Here e+ > and 
evaluating the TBA equations (6.3) at f = oo one finds that e°° = —B, and thus B = 0. 



The zero temperature limit of the equations for densities is discussed in appendix A. 2 see (A. 2. 20 1 
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Therefore the boundary is the point {^i^B) = (/iOjO) where /iq is the solution to the 
equation ei|^,^(0) = 0, with ei|„^ determined by the integral equation 



Using the identities s + K2-k s = Ki and Si\yu, = —£y® Ki this equation can be inverted 



{£y®s){Q) 



At\ 1 



dt 



i_ 

t2 



\t\>l 



(6.6) 



Evaluating this integral one finds that fiQ ranges from yUg = at weak coupling to 
lie = 2 log 2 — 2 ?a —0.6137 at strong coupling. Next for the case Qy = U , 
examining the asymptotics of the TBA equations as above yields that the boundary is 
the point {fi,B) = (0,2). Thus phases II, III, IV and V all meet here. Furthermore 
it follows that at this point ei\vu]{v) = for all v and hence there are neither particles 
nor holes of l\vw strings in the ground state here. To fill in the picture the line joining 
(/io,0) and (0,2) must be determined numerically from the coupled integral equations 

= ^ ~ B — ti\yw * Ki, ei\yw = £i\vw — 2/i + ®Qy Ki — ei|^^ -k K2 , 

parametrised by the interval Qy C X^UX^. A practical way to achieve this is to introduce 
the derivative with respect to v of these equations 

^'y ^ ^'y ~ ^'ilftf -^1' ^ ^i\vw + 4 ®Qy ~ ^il^w * -^2, 



where here we have used ^KM{t — v) = —^KM{t — v), integrated by parts, and used 



at' 



that ey{v) = on the boundary of Qy. Then for a given interval Qy the corresponding 



point {jj, B) is found by identifying the asymptotics 



-2-ii + B with 



-,t>o 



dt e'^{t) 



dte__{t) 



-,t>o 



-l\vw 



7 TBA for the extended models 



The TBA equations for the weak coupling limit of the A-model and the strong coupling 



limit of the B-model can be derived from the TBA equations (4.20 4.22) as discussed in 



sections 5.1, 5.2, Now we wish to investigate the integrable extensions of these models, 
presented in section 2.4 To this end we review the necessary modifications to the 



analysis of the previous four sections. The notation of writing in script kernels which do 



not depend on u due to a rescaling of their argument is used, see (5.1). 
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7.1 Bethe and TBA equations 

Let us first consider the models without the extension, so g = 0. To handle the Bethe 



equations (3.3) correctly it is necessary to rescale v — )■ v/u, w — )■ w/u and then take the 



appropriate limits of u. Then they become^ 



M 



a 



N+ N. M 

nWk - z+j + t YjWk- z_j +t Yj Wk - Wi + 2t 
:|| := M ^, k = l,...,M. 

M . 

e^p-.^^ = (-1)^ = rr ^-'^ "^^'^^ A; = 1, . . . , iV. 

J- J- r , — 11) ■ — 1 



. Z^^k -Wj-i 



5 ' — 5 



where a = —1 for the extended A and opposite A models, and a = 1 for the extended B 
and opposite B models. The dispersion relations for the extended A- and B-models are 

S^^[v) = -2cosp+{v) - 2 = -AnJCiiv), £^''{v) = -2 cosp_(i;) - 2 = , 

£Boo ^ _2 cos p+ (f ) = A-KlCiiv) - 2 , S^°^ (?;) = -2 cos j9_ (f ) = 2 , 

and the respective dispersion relations for the opposite models have the opposite sign. 
In the limiting process the ?/-roots decouple into two sets 

: |?/_|_ = x{uv) : Xm{v) > 0} U {y_ = — ^ y : Xm{v) < 0} , 

p_ : {|/_ = — r : Xm{v) > 0} U {?/+ = x{viv) : Xm{v) < 0} , 

with the momenta as presented in the Bethe equations above. The corresponding f-roots 
are labelled 2+ and Z-. Note that the z_-roots have been frozen, while the 2;+-roots retain 
non-trivial momentum dependence. As the Hubbard term commutes with the respective 
Hamiltonians, the number of doubly occupied sites is conserved and should be either 
or iV_. It is clear that it must be A^_, the number of frozen roots. Indeed this is 
understood as 2;+-roots being a first level of electrons that occupy empty sites, and that 
z_-roots are an extra level that combine with a 2;+-root and a w-root to make a site 
doubly occupied. As a consistency check note that there are at most M z_-roots as each 
must satisfy the same polynomial of degree M, and as 2M < + , that for each z_ 
root there exists a corresponding z+-root and w-root. 

Next we make the observation that M|fw-strings are composed of M + 1 2;+-roots 
and M — 1 2;_-roots. Let us illustrate this for the case of the l|fw-string. Recalling 
equation (4.11), yi = x{v + i) and y^i = , it is clear that the two |/-roots are of p+ 



type. For longer strings the logic extends naturally. This has an important implication. 



^^Thesc Bethe equations coincide with those derived in [37] with the BFFB grading. 



32 



The l|vw-string here corresponds to a bound pair of a spin- up electron and a spin-down 
electron, each on a singly occupied site. Indeed, the Hubbard interaction favours singly 
occupied sites and thus extending the models can favour electron pairing. 

Now consider how the coupling constant g enters the analysis. Since z+ corresponds 
to a singly occupied site the Hubbard term shifts its dispersion relation as S+ ^ £+ — 2g. 
Similarly the dispersion relation for Z- is shifted as ^ S- + 2q because it corresponds 
to a doubly occupied site. Finally since an M|f ly-string is composed of M + 1 2;+-roots 
and M — 1 2;_-roots its dispersion relation is shifted as Sm\vw £m\vw — 40. 

Despite apparent differences between the extended A- and B-models, the Bethe equa- 
tions above allow an exact equivalence to be seen. Note that they are identical for lattices 
of even length L. Moreover shifting g — t- g — 1 takes the dispersion relations for the B- 
model to the dispersion relations for the opposite A-model. Similarly, shifting g — )■ g + 1 
takes the dispersion relations for the opposite B-model to the dispersion relations for 
the A-model. Thus the spectra of the related models are identical for lattices of even 



length. At the level of the Hamiltonians (2.21 )-(2.24), this implies the following unitary 



equivalences as discussed in section 2.4 



H^°(g) ~ H^-(0 - 1) , H^o(fl) ~ H^°°(0 + 1) . (7.1) 

In particular it follows that the EKS model and the extended B-model are the same up 
to a shift of the coupling g by 1. 

Now let us discuss the TBA equations. The eigenvalue of the free energy, the analogue 
of equation ( |A.2.21 ), is 



/(/i, 5, T) = j Av{£+ - /i - 5 - 20)p+ + j Av{£_ -fi-B + 2q)p_ 

CO „ oo „ 



(7.2) 



M=l M=l 

with the dispersion relations for M|f ly-strings given by 

and the signs are changed for the opposite models. The effect of g is to shift the chemical 
potential, differently for different string types. The canonical TBA equations are thus 

1 -I- 1 

log3^+ = + log — — — -k /Cat, 



£^-fi-B + 2Q 1 + 5^ 

— h log 

T ^ 1 + 



log y. = ' + log ^ /c^. 



1 ^; £M\vw-2Mp-Ag ^ \ ir ii^^^r 
^ogyuivw = h log 1 + — ) -k ICnm + log — — Y~ * 

log yM\w = h log 1 + - — -k ICnm + log — — f -k iCm ■ 

1 V J^^fl7,>/ 1 + ^ 



(7.3) 



3^- 
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Ao 


Ao 


Boo 


Boo 






-47r/Ci 




47r/Ci - 2 


2 


-47r/Ci 












2 




-2 






-47r/C2 


47r/C2 


47rX:2 - 4 


4 


-47r/C2 




E'\\vw -5 


-4:715 


47rs 


47rs 




-47rs 






AtxKi - 47rs 


47rs - 47r/Ci 


47rs - 47r/Ci + 4 


47r/Ci 


- 4 vrs - 4 




{J, ca) 


(2,1) 


(-2,-1) 


(-2,0) 




(2,0) 



Table 2: Driving terms of the TBA equations (7.4), and constants for the free energy. 
Simplifying these equations as in appendix A.2| they become 

log 3^4 



h log TT^ -k 3 



T 



1 + 3^iK 



logj^- 

log yM\vu 



£- - £l\vw * 5 + 40 1 + 

h log * 5 . 



T 



1 + yl\^ 



1 + 



ImN log (l + yN\vw) *5 + 5 Ml log ^ ^ y -kS 



(7.4) 



lo, 



€ yM\w = hlN log (1 + yN\w) * S + 5mi log J_ * 5 



with the familiar large M asymptotics 

^OgyM\vw _ 2/i 



M^oo M 



T 



log3^M|«, 2 5 
lim — — = — — 

M^oo M T 



The driving terms in these equations are simplified and presented explicitly in Table |2} 
Taking into account that 



— = -27r/Ci , -^ = 0, — — 
av av av 



-2txKm+i{v) 



one finds that the minimized free energy 7.2 



IS 



M+1 



(7.5) 



Computing the infinite sum in the second term as in appendix A. 2 one gets 
/ = -/i-20-2cA + Jlog2 

-T(log(l + ^)*5 + log(l + ^)*/C2*s + log(l + 3^i|^^)*/Ci^s) 

where J and ca are given in Table |2] and they follow from the identities 

2713 7^ /Ci = log 2 , 27rs /C3 = 1 — log 2 . 



(7.6) 



(7.7) 
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The free energy can be also written with dependence on the l|w-string 

f = -f^ + Jlog2 - T(log(l + y+)*5 + log(l + y.)*}C2*s + log(l + 3^11^) * /Ci *s) . 

Since the extended opposite A-model is the EKS model the TBA equations and free 
energy above describe its equilibrium states in the thermodynamic limit. Our canonical 
TBA equations look very different from those derived in [2H], because we used the 
BFFB grading for the Bethe equations while [28] used the BBFF grading. The sets of 
equations nevertheless are equivalent. Identifying the functions a , /3 , 7 used there with 
our Y-functions as 

Oil = i/y+, = y-, aM+i = yM\w , A = i/3^ii™ > im-i = yM\vw , (7.8) 

we find that our simplified equations for 3^+, yM\w , M > 1 and yM\vw , M > 2 agree with 
the equations (3.13) from [2H], and our equation for y^/y_ agrees with their eq.(3.10) 
combined with eq.(3.8) for n = 1. Finally the canonical equation for 3^i|„i„ follows from 
the combination Ki *log(l + ai) + {6 + K2) *log/3i and their eqs.(3.8, 3.9),^^ 

Let us also mention that the simplified equations for and yi\vw are presented in 
[29] in eq.(5.7), and we match our equation for 3^i|^^ with theirs. The equation for 
however in [22] does not contain the crucial coupling dependent term 40/T. It seems 
that to get the term from the quantum transfer matrix approach one should specify the 



large f -behaviour of 3^+/3^_ which is determined in appendix A. 4 



7.2 Strong coupling limit 

There are two natural ways to take the large q limit of the extended model TBA equations 



(7.3). Taking the strict — )■ cxd limit with all other parameters fixed will lead to a half- 
filled model, due to the dominance of the Hubbard term favouring singly occupied sites. 

Alternatively one can take /x — ?■ — 00 simultaneously with g — )■ 00 and retain the 
possibility for empty sites in the Hilbert space. In particular, if one redefines the chemical 
potential as jl through 

/i = /i + 2fl, (7.9) 

the effect is that the models have been extended with the term 4^ X]j=i ^'.t^i.^' rather 
than the term AqV^. Then the extension does not give a cost to empty sites. Let us 
proceed with the discussion of this limit with the redefinition of the chemical potential 
and return afterwards to the strict g — )■ cxd limit. Note then that — )■ 00, and 



yM\vw — )■ C)0 for M > 2, as — 7- 00 with /i fixed. The equations (7.3) thus take the 



^■^We thank Fabian Essler for matching the equations for y+jy^ and 
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following simplified form 



log 3^4 



s+ — i^ii^u) ^5 1 + yuvw 

h log * 5 , 



T 



+ log 1 



1 + yi[. 
1 



^ /C2 + log 1 



log3^Af|«; = hiN log (1 + yN\w) *5 + 6mi log (l + :;;7- 



free energy (7.5) simplifies to 



which are the TBA equations of the supersymmetric t-J model [30]. The expression for 

^/Ci- Tlog( 1 + ^^ ) ^/C2. (7.10) 



/ = -Tlog 1 



1 



3^4 



yl\^ 



Now consider the g — )■ 00 limit with fx fixed. From the TBA equations (7.3) and 
(7.4) note that — )■ 00 and yi\vw 0, and thus the simplified equations become 



log 3^4- = 'y'""^^ -iog(i + yi|.)^5 

log yM\w = hiN log (1 + yN\w) * 5 + 5 Ml log (l + :;;7 

^ J^4 



where yM\vw have decoupled. These can be simplified further by relabelling: = 1/3^4 
and yM = yM -i\w for M > 2. Then the Y-functions form one set 



log yu = 
and the free energy is given by 



T 1 \ ^ 2n J s 

ImN log (1 + yN) * 5 + dMl rj. , 



f = - ^l-2Q-2cA + J\og2 - Tlog(l + 3^1) *s 



(7.11) 



This is none other than the set of TBA equations for the XXX-spin chain (5.8). Com- 
paring the values of J from Table |2] it follows that the g 00 limit of the A-model 
and opposite B-model is a ferromagnetic spin chain, while the B-model and opposite 
A-model become the antiferromagnetic spin chain. 



7.3 Zero temperature phase diagram 

Here the zero temperature phase diagram of the extended B-model is obtained for the full 
range of coupling G M, see Figure [1} At the end we will comment how the corresponding 
phase diagrams for the other three extended models can be obtained from it. 

Let us break the analysis into two parts. First consider the case g > 0. As in section 
[6] there are no M\vw- or M|w-strings for M > 2, and as 3^+ < 3^_ there are no w- 
particles at zero temperature either. Thus neither are there z.-particles, as can be seen 
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rigorously by examining the equations for zero temperature root densities (A. 2. 20). The 
TBA equations for the dressed energies are then 



which differ from equations (6.3) only by a shift of fi by 20. Thus the phase diagram 
will be that of the strong coupling limit of the B-model with /x shifted as /i — )■ /i — 20. 
The question that remains to be addressed is how phases III and V are separated for 
the range /i G (—20, 0). Recall that the half-filled phase is determined by the condition 
that ei\yw{v) < for all f G M, and the transition to the half filled spin polarised phase 
is when e+(f) < and ei\yu,{v) = for all f G M. From the large v asymptotics of (7.12) 



it follows that the phase boundary is the line B = fi + 2 + 2q. 

Now consider < 0. Here one can no longer say that y+{v) < 3^-(f ) for all v and so 
there is the possibility of w-particles and ^--particles entering the ground state. Recall 
however that there can be no ly-particles, and consequently no z_-particles, unless there 
are z+-particles. Note though that e+ < puts restrictions on fi and B that result in 



€i\w > and e_ > 0. Thus equations (7.12) govern the dressed energies for < also, 



and again the phase diagram is that of the strong coupling limit of the B-model shifted 
by /i — )■ — 20. Let us remark that phase V will disappear when = — 1 + log 2 and 
that phase IV will disappear when = — 1. 

The phase diagram for the EKS model (extended opposite A-model) is obtained from 



the one for the extended B-model by shifting the coupling above by 1, see (7.1). For 
B = our results agree with . The corresponding phase diagrams for the extended A- 
and opposite B-models are obtained by using the relationships between opposite models 



discussed in appendix |A.3 



8 Conclusion 

Given the wealth of work that has been done on the Hubbard model there are many 
possible directions in which one can take the A- and B-models. To understand better 
the physical properties of the models it is desirable to analyse the spectrum of excited 
states and the finite-size corrections at zero temperature. Taking into account the rather 
involved form of the Hamiltonians it would be interesting to consider the continuum limit 
of the models. 

The sets of TBA equations we obtain depend on infinitely many functions. It would 
be useful to apply the fusion hierarchy approach to the thermodynamic limit and derive 
a finite set of nonlinear integral equations from the quantum transfer matrix [HT| H] • 

In this paper we have restricted our focus to the parity-invariant Hubbard-Shastry 
models, but it might be of interest to analyse the general parity-breaking models. Such 
chiral models could exhibit interesting behaviour, see for example [33] . 

It is worth noting that there are other integrable extensions of the Hubbard model, 
in particular the Alcaraz-Bariev set of models These correspond to a quantum 
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deformation of Shastry's R-matrix [M], and so provide anisotropic deformations of the 
symmetries of the A- and B-models. This might allow one to introduce spin-orbit cou- 
pling to the models. 

The nonintegrable extended B-model interpolates between the Hubbard, EKS and 
t-J models. Thus it might be a useful toy model for studying transitions between the 
various phases they exhibit. In particular the t-J model is the established starting point 
for the study of the high temperature superconductivity [35], |36j, while the Hubbard 
and EKS models are also mentioned in this context. Given the high symmetry of the 
extended B-model it might capture interesting physics that its t-J limit does not. 

Another direction is to construct models for the Kondo lattice [SZl ES]- These are 
models with eight states at each site, where the pseudo-vacuum carries spin-1/2 and 
so the site vacuum can be labelled | t); I i)- The spin of the itinerant electrons can 
couple to the spin at the site and the models are used for the study of heavy fermion 
systems. Such models arise naturally from the R-matrix for the 8-dimensional atypical 
representation of su(2|2)c constructed in 

Finally it is of interest to investigate how the models behave in higher dimensions. 
Our analysis has been largely restricted to one dimension as it is only there that the 
models are integrable. Higher dimensions however are more physically relevant and may 
also exhibit rich behaviour. Given that the models possess high symmetries there is also 
the possibility that some exact statements can be made. 
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A Appendices 

A.l Conventions, definitions and notations 

Conventions for obtaining spin chain Hamiltonians from R-matrices 

To compute the Hamiltonians of the models we have used the su(2|2)c-invariant R-matrix 
in the form given in [39j. Up to an overall scalar factor it is related to Shastry's R-matrix 
by a similarity transformation. The R-matrix depends on two independent parameters, 
and it can be thought of either as a function ^{xf^x"^) of the four parameters xf ,xf 
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subject to the two constraints 



1 



X I 



Aiu, k = l,2, 



X,, 



(A.1.1) 



or as a function R(-2i, Z2) where Zi , Z2 parametrize the tori defined by (A.1.1 ). One can 
show that xf ,xf and Il{zi,Z2) are meromorphic functions of Zi and Z2 [39j. For this 
reason it is safer to use the z-parametrization for computations. The Hamiltonian is 
given by H = ^^^^ where Hjj+i is computed as 

H^2 = ^^fn{z)P^2^^R,2{z^,z)l^=,, (A.1.2) 

where the normahzation constant Nuiz) is chosen so that the dispersion relation would 
be given by (3.6), and since the spin chain is homogeneous it is sufficient to compute 
only H12 . We find the following normalization constants 

(H) z = uji/2 - UJ2/2 , Afuii{z) = 2u, 

(A) z = -UJ2/2 , A^ha(^) = 2uVl + , 

(B) z = uJi/2, A/hb(^) = 2u/Vl +u2. 

Here coi = 2K{—l/u^) and UJ2 = 2iK(l + l/u^) — 2K{—l/u'^) are the real and imaginary 
half-periods of the torus defined by (A.1.1). 



One can also obtain the same Hamiltonians as 

d 



H 



12 



u(x 



12 



Ri2{xf,x^^ 



dxi 



where Xi is considered as a function of Xi- One then gets 

x~^ , , 



^fuix^) 



(X 



X 



x~ 



1 — X ^ 
1 - X+2 ' 



(A.1.3) 



(A.1.4) 



and therefore 
(H) 

(A) x-^ 

(B) x^ 



l/x" = 00, 

-1/x- = i{u + Vl + u2), 
—X" = i(u + vTTu^), 



A^HH(a;=^) x+/x~ , 
^fiiA{x^) = 2uVl + u2 



HB(X 



2^/T 



The resulting H12 is represented either as a 16 x 16 matrix or as a differential oper- 
ator and can be rewritten in terms of the fermion creation and annihilation operators 
c|, ^ , Cfc^f- ,cl^, Ck^i by using the identification 



E:^^ — Eo 



Ea 



= -^31 — -C/24 , = -[^13 — -^i2 , = -C/41 + -E23 , Ci = Eii + E32- (A. 1.5) 

c| (g) /4 O cj_^ , C| (g) /4 O Ci^^ , c| (g) /4 O c\^^ , C4. (g) /4 O Ci_^ . 
P-h®c\-P O P-hmyP ^ C2,t, /^■/4®c|-/^ ^ cj^^, P-h®CyP ^ €2,4. 

Here are 4x4 matrix unities, /4 is a 4 x 4 identity matrix, and P = Y2i ji^^Y"'^^ En (g) 
Ejj is the graded identity where ei = e2 = 0, es = €4 = 1. 

Finally, to get the kinetic term in the canonical form (2.2) one should perform the 
canonical transformation (2.13) with a = —7i/2 for the Hubbard model, a = ir for the 
A-model, and a = —it/2 for the B-model. 
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Matching the notations and conventions 

Most of our notations and conventions come from [IQ], and here we compare them to 
those of [4J. 

In the Bethe ansatz we denote particles momenta as pj and auxiliary roots as Wj, so 
they are related to the ones in [1] as pj -f-)- kj, wj f-)- Aj. 

In the string hypothesis a M|w-string is a A string of length M, a ty-particle is a 
A-string of length 1, a M|fw-string is a k-A string of length M, and ^/-particles could 
have been called fc-particles. 

In the TBA equations the Y-functions are related to the ones in [4j as Ym\w ^ Vm, 
Ym\vu, ^ v'm, y-(sin(A;)) o ((k) , \k\ < tt/2 and Y+{sm{k)) o ((k) , \k\ > 7r/2. 



Kernels and S-matrices 



(4.16 


4.18 


) and ( 


A.2.18 



equations involve convolutions with three types of kernels which we list below 



s{v) 
Km{v) 
Kmn{v) 



2Tiidv^^ ^""^ 4ucosh'^^' 



1 d 



27ci dv 
1 d 



log Sm{v) 



2u 

uM 



, , ,TTV in. 

b[v) = — tanh 



Sm{v) 



V — iuM 



27ci dv 
Smn{v) = Sm+n{ 



Note that the kernels satisfy the following identities 

1 ★ s = - , l-k Km = 1 , Km * Kj^j = Km+n , 
Ki - s-k K2 = s , Km+i - s-k Km - s-k Km+2 = 



(A.1.6) 
(A.1.7) 



v + iuM 

M-l 

log Smn{v) = Km+n{v) + Kn-m{v) + 2 ^ KM-M+2j{v) , (A. 1.8) 

M-l 

v)Sn-m{v) W SN-M+2j{vf = Snm{v) ■ 



(A.1.9) 



(A.1.10) 



A. 2 Deriving the TBA equations 

String hypothesis 



The string hypothesis for the models has been formulated in section |4.1[ and according 
to it in the thermodynamic limit every solution of the Bethe equations (3.3) corresponds 
to a particular configuration of Bethe strings and consists of 

1. ?/+-particles with rapidities v+^k 

2. ?/_ -particles with rapidities V-^k 

3. Nm\vw M|t>w-strings, M = 1, 2, . . . , 00 with rapidities VM,k 
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4. Nm\w M|w-strmgs, M = 1, 2, . . . , oo with rapidities WM,k- 

The numbers N^, Nm\vw, Nm\w are called the occupation numbers of the root 
configuration under consideration and they are related to A^, M as follows 



N = N+ + N^ + Y.2QNq\,^, 



g=i Q=i 



Bethe equations for string configurations 

Next we rewrite the Bethe equations in terms of the string configurations. The S- 
matrices presented in equations ( A.1.7[[A".1.9 ) are used. The Bethe equation for a y± 
particle with rapidity v±^k is 



oo ^M\vw 



N=l j=l 



M=l i=l 



where the p± are given by equations (3.2) and (4.8). The Bethe equation for an M\vw- 
string with root VM,k is 



oo ^N\vw 



N+ Af_ 

= e^PMi„.K.,.)L -Q SM{vM,k - v+,) n SM{vM,k - v.,j) n n SMN{vM,k ' ^^-12) 

i=l j=l 



N=l 1=1 



and the momentum Pm\vw is given by (4.12). Finally the Bethe equation for an M\w- 
string with root WM,k is 



oo ^N\w 



n n ^MN{wM^k - Wn,i) = Y\_ SM{wM,k - W Sivl{wM,k - W-,j)(A.2.13) 



N=l 1=1 



1=1 



TBA equations for densities 

We introduce densities p{y) of particles, and densities p{v) of holes in the standard way 
and obtain the following TBA equations for the densities by taking the log derivatives 
of the Bethe equations for the string configurations 



P+ + P+ 

P- +P- 

PM\vw + PM\vw 

Pm\w + Pm\ w 



1 dp+ , . 

— 27r ~(hj ^-^ ^^^^l^"" ~'~ P^\^ ) 

M=l 

1 dp^ 1^ , s 

2n~dv~ ~'~ ^ \Pm\vw + Pm\w) , 



(A.2.14) 
(A.2.15) 



PN\vw KM^ip+ + P-), (A.2.16) 



N=l 



M=l 
1 dpM\vw 

2tt dv 

oo 

Pn\w 

N=l 
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(A.2.17) 



Here we have taken into account that for all the models ^ < and > 0, and 

dv dv ' 

"^^^^'"^ < 0. The convolution i is defined as 

dv 

g-kh = g-kh for the Hubbard and A — models , (A. 2. 18) 

gk:h = g-kh for the B — model , 



and the kernels are defined by (A. 1.7) and (A. 1.8). It is worth noting that 



l*(p+ + p_ + p+ + p_) = 1. (A.2.19) 

Let us briefiy consider the zero temperature limit of the equations for densities for 
A- and B-models. As is discussed in section |6] in this limit Ym\vw > 1 for M > 2 and 
yM\w > 1 for M > 1, and therefore the corresponding densities of particles go to zero. 
For the remaining string types there are intervals Qk where the densities of holes vanish, 
and outside which the densities of particles vanish. Then the equations above simplify 

1 dj)—. 

P- =^^ + ^i*Qii.„Pi|.«;, v^Q_, (A.2.20) 
1 dpi\vu, J. ^ 

Canonical TBA equations 

The free energy in terms of the string densities takes the following form 



/(p, B,T)=J dvi£+ - p - B)p+ + J dv{£. - p - B)p. 

(A.2.21) 

+ 5Z / df(£^M|t,«, - 2Mp)pM|,;«, + ^ / dv2M B pm\w -T s. 

M=l M=l 

Then the entropy per site is given by 

s = dv is (p+) + s (p_) + ^ {s{pm\vw) + s{pm\w)) I , 

\ M=l / 

where s(p) denotes the entropy function of densities of particles and holes 

s(p) = plog (l + ^) + Plog (l + ^) • (A.2.22) 

The TBA equations are obtained by minimizing the free energy and introducing the 
Y-f unctions 

^ _ P± ^ _ PM\vw ^ _ pM\w ( \ n 

i± — , 1 M\vw — 5 1 M\w — 5 (^A.Z.Zdj 

P± pM\vw pM\w 

42 



they can be written in the form 

log y± = + log ^ K^^ (A.2.24) 



Yj 



Nliv 



logFMi.. = ^^^1--^^^ + log (l + -LJ] ^ K^^^ _ log (i + ^)®Km, 

(A.2.25) 

log Ym\u, = + log (^1 + * i^jvA/ - log (^1 + ^^®Km, (A.2.26) 

where we sum over = 1, 2, . . . and recall that Y±{v) = Y{v±iO). We refer to the TBA 
equations written in this form as the canonical ones. The minimized free energy is 

The infinite sum in the second term can be simplified as will be discussed shortly. 
Simplifying the TBA equations 

The above equations can be rewritten in a simpler "local" form. To achieve this we 
introduce the following kernel 

{K + (v) = 5MN5iv) - Imn siv) (A.2.28) 

with Imn = Sm+i,n + 5m-i,n , M >2 and Iin = that is inverse to Knq + 6nq- 

00 00 

J] (K + 6)-J^ * {Knq + 6nq) = i^Q^ + ^Qn) * + ^yNM = ^MQ , (A.2.29) 

N=l N=l 

which can be also written in the following convenient form 
00 00 

+ 6)-;^ ^K^Q = J2 Kqn ^{K + 6)-\, = Imq s . (A.2.30) 

N=l N=l 

Here and in what follows we often use the convention SMNi^j) = ^mn^{v) when this will 
not cause confusion. This obeys the useful identities 

00 00 

^ (ir + ^Kn = Y,Kn^{K + 5)-^\, = 6mi s , (A.2.31) 

Af=l Ar=l 

Sn\vw * (i^ + 5)^M = -^iM ^®s = Sim {£+- £-)^s. (A.2.32) 
Now we can proceed with simplifying the canonical TBA equations. Subtracting the 



TBA equations (A.2.24) for y-t-particles gives 



log^ = (A.2.33) 
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Then applying the inverse kernel to the equations for the vw- and w-strings and using 
the identities above and (A. 2. 33) gives the simplified TBA equations: 

log yiVl\vw = hlN log (1 + Fiv|^«,) * s - 6mi log {1 + Y)®s, (A.2.34) 

log Ym\w = hiN log (l + Y]^\^) -ks- 6mi log (^^ + ® , (A.2.35) 

where we took into account that {K + 6)~^ annihilates the fi and B dependent terms. 

The simplified TBA equations can be used to compute the infinite sums appearing 
in the canonical equations for ^/-particles. We show how to do this on the example of 
vw-strings. Since a similar sum appears in expression (A. 2. 27) for the free energy we do 
the computation by assuming that we are given a kernel Jf^f satisfying the identity ^ 

oo 



(A.2.36) 



N=l 



where 5J(f is any kernel. We want to compute the following sum 

1 



M=l ^ ^M\vwJ 



M 



(A.2.37) 



To this end (A.2.34) is rewriten in the form 

hlN log (l + 



^Og M\vw 



Imn iogF;;^^*s 



N\vw 



^ 1 
)*S-(5Ml(^ + log (1 + ^) 



where on the l.h.s. of this equation we replaced Ym\vw by Y^^^,^,^ = Ym\vw^'^^''^^'^ to make 
sure it asymptotes to 1 at large M. This does not change the equation because the 
kernel s integrates to 1/2, that is 1 * s = 1/2. This replacement is necessary because 
at the next step we multiply this equation by J^m-, take the sum over M and use the 
identities (|A.2.36l) and (|A.2.32l) to get 



\ogYi;i^5Jf = log 1 + 



Y, 



M\vw 



+ 



* JTm - log I 1 + 

1 



1 



T 



log 1 + 



Y 



Yi\, 

S -k .J^i . 



i.5Jf (A.2.38) 



From this equation and the definition of Y^^^^^ we immediately obtain 



log 1 + 



Y, 



M\vw 



-k J(fM = log (l + Yi\^^) -k 5Jf 



^1 \vw 2^ 

T 



+ log 1 



Y 



s^JTi. (A.2.39) 



This method was apphed in |40j to more general kernels but under an assumption of Y^f-functions 
approaching 1 at large M. 
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A similar formula can be derived for the infinite sum with YMiw-functions 



log ( 1 + ) * JTm = log (1 + * 5Jf 

V iM\wJ 



2B 



* 5,J^ + log 1 



Y 



S -k . 



In our case ,J^m = Km, = s, and subtracting these two equations one obtains 

, ^ , 1 + Y\\vw , fJ' + B Sil^u, 

log 1 — -k Km = log , . — s -\ ic s . 



1 + 



Yn.. 



1 + 



The infinite sum is substituted for in (A. 2. 24) to give the simplified equation for V+Yl 
log Y^Y_ = — h 2 log TT — * s . 



T 



1 + 



Combining the equations for Yj^^/Y^ and Y^Y_ gives 

logF = 

Simplifying the free energy 



h log -k S 



(A.2.40) 



Eq.( A.2.39 ) can be also used to compute the sum appearing in the free energy expression 
( |A.2.27D . First we rewrite ( |A.2.39D as 

log ( 1 + ^T—] ^J^M = log (1 + * + ^ * + log (1 + F) ® s * JTi . 

\ yM\vw J j- 

(A.2.41) 

Then we notice that 

'dp+ dp_ 



and therefore in this case J^^a 



M 



r dp _ 

df 



dv 



dv 



dv ' 



SJ^ = -s®^. Thus eq.( |A.2.4lD becomes 



log 1 



Y, 



M\vw 



dpM\t 

dv 



1 A \ dp 2u 

log (1 + *s® — - — vts 



dp 
df 



+ log (1 + r) ® s * 



dv 



Taking into account that 1 t^t s t^t — = — vr , one gets 



/ = -A.-£log(l + p 



dp T ^ . dpiu 

+ — log 1 + r ® s * 
df 27r dv 



(A.2.42) 



(A.2.43) 



^log(l + Fi|„^)*s®^ 
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Finally using eq.(A.2.40) and identity (A. 3. 50) the free energy can be rewritten as 



f 



dp 
dv 



T f 1 



/dp 
\dv 



S -k 



dv 



T d p 

— log(l + Fi|,^)vcs® — 



or with dependence on the l|iu-string in the from (4.15) 



(A.2.44) 



A. 3 Attractive vs. repulsive models 

For any given model one can consider H or H as the Hamiltonian. Here it is useful to 
distinguish the two cases by the sign of the coefficient of the Ylj ^jt^ji term, calling the 
model repulsive if the coefficient is positive and attractive if the coefficient is negative. 
Throughout the text the repulsive case was singled out as it appears to be the more 
physically relevant. A primary application of the models presented is to systems of spin- 
I electrons and these are driven by the Coulomb repulsion. Nevertheless, it may also be 
of interest to examine the attractive cases. 

The TEA equations for the attractive Hubbard model were derived in [H]. Here we 
show how the TEA equations and free energy for the attractive Hubbard-Shastry models 
can be straightforwardly obtained from those for the repulsive models. Up to a unitary 



transformation, see (2.15), an attractive model differs from the repulsive one by the sign 



of the Hamiltonian. They share the same set of Eethe equations but have opposite signs 
of the dispersion relations. Thus the only change to the thermodynamic Eethe ansatz 
analysis is the sign of the dispersion relation and therefore the TEA equations for the 
attractive models are 



logF 



- + log 7^ ^ S . 



(A.3.45) 



loffF/ 



M\w[ 



T " 1 + 

Imn log (l + * s - 6mi log(l + Y)®s , (A.3.46) 

Imn log (l + * s - 6mi log(l + l)®s , (A.3.47) 



and they are supplemented with the large M asymptotics 



lim 

Af-s>oo 



log Ymk 



2j2 



lim 

M-s>oo 



log Ym\w _2B 



(A.3.48) 



M T ' M^oo M T 

Here and in what follows tilded quantities refer to the attractive models, and untilded 
ones to the corresponding repulsive ones. The free energy of the attractive models 
are given by either (A.2.44) or (4.15) with the replacement £ — > —S, Si\yu} — > —Si\vw, 



Y -^Y, Yi|^ Yii^ and Y^^^ Yi|^^. 

The TEA equations are identical to the ones for the repulsive models (4.20 - 4.22) 
under the identification 
1 



Y 



Y'' 



YmIvw — Ym\w , yM\w — yM\vw , fJ' — —B , B — —fi . (A. 3. 49) 
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By using this identification one can relate the free energy of the attractive models to 
those of the repulsive one. Taking into account the identities 



s 



£l\vw * S® 



dp 



where £.„/2 is an electron energy at p 
£ 



dv 

7r/2, that is £t,/2 



dp 
dv 



£-k/2 



(A.3.50) 



-2u for the Hubbard model, 



-it/2 



2-\/l + V? for the A-model and £^^/2 = for the B-model, one finds 



f{fx,B,T)=L 



£ 



n/2 



(A.3.51) 



Since £^7^/2 depends only on u this will not affect the physical properties as they are 
related to derivatives of /. In fact £tj/2 can be removed by redefining the Hubbard 



interaction potential (2.3) as follows 



n 



1 



:)(nfc„ 



:) 



(A.3.52) 



which is the potential used in [4j. 

The same relation between the free energies of the attractive and repulsive models 
can be also found by using the Woynarovich transformation [12] which is a superposition 
of the parity transformation and the partial particle-hole transformation of the spin- 



up electrons, see |1] for detail. Under the transformation the Hubbard potential (2.3) 



transforms as Vj^^ — )■ — Vj ^fc — | , and all the other kinetic and potential terms just change 



their signs. The relation (A.3.51) then follows from the definition (4.1) for free energy 



and that N — )■ — 2S^ + 1 and — )■ — |N + ^ under the Woynarovich transformation 
A. 4 Algebraic limit of TBA equations 



As was discussed in section 5.1, in the weak coupling limit and for values of v on 
and the TBA equations of the Hubbard and B-models reduce to the set of algebraic 



equations (|5.4||5.6|). The general solution to these equations is well-known, see e.g. 

sinh(/„ + M) 



and can be written in the form 

'sinh(/„^ + M)c 



1 + Ym\vw 

1 + Y- 



sinh c 



/ sinh fy^c 
I sinh 



1 + Ym\w 



1 + 



sinh Cy 

sinh fy,Cu 
sinh c, 



Here the constants Cy^ and are fixed by the large M asymptotics (|4.24|) of Y-functions 

B 



c, 



T 



/3c 



equations (5.4) 



sinh fy^(3c 
sinh fy,l3s 



Cy, = = f3s , and the two remaining functions are then found from 

e^^ sinh /3c (cosh /3c + e~^^ cosh Ps) 
\/ 2 cosh /3c cosh cosh /3a + cosh^ /3c + cosh^ /3s + sinh^ /3a 

sinh /3s (cosh /3c + e~^^ cosh /3s) 
■\/ 2 cosh /3c cosh /3s cosh /3a -I- cosh^ /3c + cosh^ /3s + sinh^ /3a 
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with /5a = ^- The exphcit solution to the above TBA equations is then 



1 + ^M\vw 



1 + Ym\w 



e^/^'^csch^/Sc (2e-^-^ cosh sinh(M^c + Pc) + e'^^^^ sinhM/S^ + smh(M^c + 2/3^))^ 
4 (2 cosh /3c cosh /3s cosh /3a + cosh^ /3c + cosh^ f3s + sinh^ /3a) 

e'^f^^csch^Ps {2e-^^ cosh/3c sinh(M/3, + /3s) + e'"^!^^ smh(M/3s + 2/3^) + sinhM/J^)^ 
4(2 cosh /3c cosh Ps cosh /3a + cosh^ /3c + cosh^ /3s + sinh^ /3a) 
4e-^^ cosh /3c cosh /3s + 2 cosh 2/3c + e'"^^^ + 1 



4e-/3A cosh /3c cosh /3s + e-'^P^ (2 cosh 2/3s + 1) + 1 ' 



(A.4.53) 



It is not difficult to see that the same formulae also describe the large v asymptotics 
of the Y-functions of the B-model for any u because at large v Y-functions asymptote 



to constants and therefore the simplified TBA equations (4.20 4.22) reduce to eqs.(5.4- 



5.6) with A replaced by —2. Thus the asymptotic values of Y-functions are given by 
(A.4.53) with /3a = —2/3. Furthermore these formulae determine the large v asymptotics 
of the extended A- and B-models with A = — 2g for the A- and opposite A-models, and 
A = —2 — 20 for the B-model and A = 2 — 20 for the opposite B-model. 



A. 5 t-J model from large on-site repulsion 

In this appendix we consider the effect of large on-site Coulomb repulsion on the A- 
model and extended B-model. To this end it is useful to present the Hamiltonians of 
the models in the form 

H(7, C, K) = E (t„« + 4-, V« „ + C (V« „ + V- , - , + ,) ) . 

Ti,fe = - 5^ (S>cfc,, + cl ^^Cj^a) (1 - (1 - (nj,-a - nfe,_^)^) , (A.5.54) 



where for the extended B-model 7 = 0, C = 2tanhz/ and Hj = 1/ coshi/, while for the 
A-model 7 = cosh ly, ( = —2/ cosh u and kj = (—1)-' tanh u. 

This form of the A-model Hamiltonian is obtained by applying the unitary trans- 
formation (2.30) to the Hamiltonian (2.18) with a = 7r/2. Notice that as a result 



of the transformation the charge su(2) generators take the same form as the ones for 



the Hubbard and B-models (2.4). The correlated hopping term however becomes site 
dependent. 

In the large 7 limit the Hubbard term dominates and the models are half-filled. It is 
useful to shift the chemical potential /i — )■ /i + 27 to study the less than half-filled models, 
as this will remove the term proportional to N from the Hubbard interaction ( |2.7 ). Then 



let us write the Hamiltonian (A.5.54) as H(7, (, k) = 47 ^^^i nj^^nj^^ + H(0, (, k) and 



note that in the strict 7 — )■ 00 limit any state that does not contain doubly occupied 
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sites is a ground state and there is huge degeneracy. Taking the sub-leading terms into 
account in ^ perturbation theory in a standard way, see e.g. section 2. A and particularly 
equations (2. A. 26, 2. A. 30) of the book the effective large 7 Hamiltonian is 

1 ^ 

H*-^(7, C, k) = Po H(0, C, /t) Po + ^ 5^ Po H(0, C, k) n,, tn,- ^ H(0, C, k) Pq , 

where Pq = 11^=1(1 ~ '^j.t'^j,^) projector onto the model's Hilbert space of dimen- 

sion 3^ where there are no doubly occupied sites. As happens for the similar analysis 
of the Hubbard model, the second term in the above expansion gives rise to three-site 
terms as well as two-site terms. The three site-terms however will be suppressed in the 
physically interesting regime where the model is close to half-filling and so we will ignore 
them here, as is generally done. 

The effective Hamiltonian is thus 

H-^(7, C, «:) = Po E ( - SV^.+I-^ - sVi..^.'^ + (C + (^''+1 - ^))Po , 

2 

which describes the t- J model with the coupling J = C + For the A-model expanding 
J in powers of 1/u one gets J = — 1/u, and therefore a ferromagnetic t-J model. For the 
B-model J = 2tanhz/ -|- ^ ^ , and at z/ = 00 it is the supersymmetric t-J model. 
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